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THERMAL  BUCKLING  ANALYSIS  FOR  STIFFENED 
ORTHOTROPIC  CYLINDRICAL  SHELLS* 

By  L.  K.  Chang**  and  Michael  F.  Card 
Langley  Research  Center 

SUMMARY 

A  theory  for  thermal  buckling  of  an  orthotropic,  multilayered,  stiffened  cylindrical 
shell  is  presented.  The  theory  includes  the  effects  of  eccentricity  of  layers  and  stiff¬ 
ening,  and  deformations  prior  to  buckling.  It  is  sufficiently  general  to  account  for  dis¬ 
crete  rings  and  averaged  properties  of  longitudinal  stiffening,  as  well  as  arbitrary  tem¬ 
perature  distributions  through  the  thickness  of  the  shell  and  depth  of  the  stiffeners.  Two 
computer  programs  are  described  corresponding  to  solutions  for  buckling  obtained  by 
using  finite  differences  and  determinant  plotting  or  modal  iteration. 

Computed  results  for  thermal  buckling  of  unstiffened  and  ring-stiffened  shells  are 
presented  and  are  in  reasonable  agreement  with  published  results.  The  interaction  of 
thermal  loading  and  axial  compression  in  two  large-diameter  stiffened  shells  representa¬ 
tive  of  a  launch  vehicle  interstage  and  a  preliminary  supersonic  transport  fuselage  design 
is  investigated.  Results  indicate  that  buckling  can  occur  in  both  structures  at  a  realistic 
temperature  under  thermal  loading  alone. 

INTRODUCTION 

Advanced  aerospacecraft  of  the  seventies,  such  as  space  shuttles  and  hypersonic 
aircraft,  will  probably  experience  severe  thermal  environments  during  flight.  A  renewal 
of  interest  in  thermal  stress  and  thermal  buckling  problems  is  anticipated,  since  it  is 
unlikely  that  thermal  protection  can  completely  eliminate  stresses  and  deformations 
arising  from  thermal  gradients  in  stiffening  frameworks,  bulkheads,  and  covers.  Ther¬ 
mal  buckling  problems  in  unstiffened  shells  have  been  investigated  in  references  1  to  8. 
Results  for  stiffened  and  sandwich  shells  are  given  in  reference  3  and  references  9  to  12. 
An  excellent  summary  of  existing  thermal  buckling  analyses  and  experiments  is  given  in 
reference  13.  A  review  of  these  references,  however,  indicates  that  further  studies  of 
thermal  buckling  in  stiffened,  multilayered  shells  are  needed, 

*Part  of  the  information  herein  was  presented  at  the  AIAA/ASME  11th  Structures, 
Structural  Dynamics,  and  Materials  Conference. 

**NRC-NASA  Resident  Research  Associate,  now  at  Argonne  National  Laboratory, 
Argonne,  Ill. 


Recent  advances  in  computer  technology  together  with  derivations  of  more  consis¬ 
tent  buckling  theories  (refs.  14  to  17)  permit  the  development  of  buckling  analyses  which 
accurately  account  for  thermal  stresses  and  deformations  in  stiffened  shells.  The  pur¬ 
pose  of  the  present  paper  is  to  present  an  analysis  for  thermal  buckling  of  stiffened 
cylindrical  shells,  to  describe  two  computer  programs  based  on  solutions  for  buckling, 
and  to  report  the  results  of  thermal  buckling  calculations  for  unstiffened  and  stiffened 
cylindrical  shells. 

The  analysis  was  developed  by  extending  the  theory  of  reference  17.  The  theory  is 
sufficiently  general  to  account  for  the  averaged  properties  of  longitudinal  stiffening  and 
multiple  orthotropic  layers  with  differing  coefficients  of  expansion  as  well  as  arbitrary 
temperature  gradients  through  the  thickness  of  the  shell  and  depth  of  the  stiffeners.  A 
general  outline  of  the  theory  is  given  and  buckling  solutions  suitable  for  rapid  computer 
calculations  are  presented.  Two  computer  programs  based  on  determinant  plotting  and 
modal  iteration  are  described.  Results  of  computations  for  unstiffened  and  ring-stiffened 
cylindrical  shells  under  thermal  loads  are  compared  with  results  from  existing  buckling 
theory.  Results  are  presented  to  demonstrate  the  effects  of  interactions  of  axial  com¬ 
pressive  loads  and  thermal  loads  on  buckling  of  large-diameter  cylinders  typical  of  those 
which  might  be  considered  for  an  interstage  structure  of  a  launch  vehicle  and  a  fuselage 
of  a  supersonic  transport. 


SYMBOLS 

The  units  used  for  the  physical  quantities  in  this  report  are  given  both  in  U.S. 
Customary  Units  and  in  the  International  System  of  Units  (SI).  A  table  of  conversion  fac¬ 
tors  is  given  in  appendix  A.  The  relationship  between  these  two  systems  of  units  can  be 
found  in  reference  18. 

A  cross-sectional  area  of  stiffener 

Ai,Bi,Ci  4x4  matrices  defined  after  equation  (C16) 

Cjj  extensional  stiffness  of  stiffened  shell  wall  (eqs.  (BIO)) 

Djj  bending  stiffness  of  composite  shell  wall  (eqs.  (BIO)) 

E  Young's  modulus 

Ex,Ey  Young’s  moduli  of  orthotropic  layer  of  shell 
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G 


shear  modulus 


a  length  of  cylinder 

d  stringer  spacing 

e  distance  from  reference  surface  of  cylinder  to  surface  on  which  acts 

f  function 

g  shape  factor  for  exponential  temperature  variation  (eqs.  (B5)  to  (B7)) 

h  height  of  stringer  or  ring 

i  integer 

k  total  number  of  finite-difference  intervals  along  length  of  cylinder 

I  ring  spacing 

n  number  of  full  waves  in  cylinder  buckling  pattern  in  circumferential  direction 

p  pressure  (positive  for  external  pressure) 

Pq  sea-level  atmospheric  pressure 

t  thickness  of  isotropic  cylindrical  shell 

Ae 

t  effective  wall  thickness  of  stiffened  isotropic  cylinder,  —  +  t 

u,v,w  displacements  of  stiffened  shell  at  reference  surface  in  x-,  y-,  and 

z -direction,  respectively 
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x,y,z  orthogonal  curvilinear  coordinates  (see  fig.  1) 

z  distance  from  centroid  of  stiffener  to  reference  surface  (see  fig.  1),  positive 

if  stiffener  lies  on  external  side  of  reference  surface 

A^j  structural  coefficients  defined  after  equations  (CIS) 

Ilg  strain  energy  of  composite  shell  wall  (shell  plus  stiffeners) 

potential  energy  of  loading 

a  coefficient  of  linear  thermal  expansion 

coefficients  of  linear  thermal  expansion  for  orthotropic  layer  of  shell 

/3sk>/^r»^s  temperature  distribution  factors  (see  eq.  (Cl)) 

/3i>^2v  •  -(^7  structural  coefficients  in  equations  (C4)  and  (C6) 

^sk’^r’l's  temperature  distribution  factors  (see  eq.  (Cl)) 

Vjjy  shearing  strain 

A  distance  between  adjacent  finite-difference  stations 

noo 

6(x  -  jZ)  Dirac  delta  function  defined  so  that  \  f(x)  6(x  -  jZ)  dx  =  f(jZ)  where 

^  -OO 

6(x  -  jZ)  =  0  when  x  #  jZ 

®iH  Kronecker  delta,  equal  to  zero  when  i  ^  H  and  equal  to  1  when  i  =  H 

€x,ey  normal  strains 

changes  in  curvature 

[I  Poisson's  ratio  for  isotropic  material 

Poisson's  ratios  for  extension  of  orthotropic  layer  of  shell 
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normal  stress 


|=Wa 
a 

T  shear  stress 

Subscripts: 

A  prebuckling  state 

B  small  changes  away  from  prebuckling  state  which  occur  at  buckling 

H  integer  denoting  finite-difference  station  at  which  ring  is  located 

cr  at  buckling 

i  index  of  finite-difference  station 

max  maximum 

r  rings  (circumferential  stiffening) 

s  stringers  (longitudinal  stiffening) 

sk  skin 

Superscripts: 

i  layer  index 

t  total  strain  in  composite  shell 

A  subscript  preceded  by  a  comma  indicates  partial  differentiation  with  respect  to 
the  subscript. 

Primes  denote  total  derivatives  with  respect  to  x. 
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BUCKLING  THEORY 


Assumptions 

The  theory  employed  herein  is  an  extension  of  the  theory  of  reference  17.  The 
cylinder  considered  (see  fig.  1)  is  composed  of  a  multilayered  orthotropic  shell,  stiffened 
by  uniform,  equally  spaced  rings  and  stringers,  all  having  linearly  elastic  properties. 

The  elastic  constants  of  the  multilayered  shell  are  taken  as  those  given  in  reference  19. 
Individual  layers  of  the  shell  are  assumed  to  be  such  that  their  principal  axes  of  orthot- 
ropy  are  alined  with  the  principal  directions  of  loading.  The  stringers  are  assumed  to  be 
closely  spaced  so  that  their  elastic  properties  may  be  averaged  over  the  stringer  spacing. 
The  effects  of  the  rings  are  not  averaged  but  instead  the  rings  are  considered  to  be  located 
along  the  length  of  the  shell  at  uniform  spacings.  The  effects  of  the  individual  (discrete) 
rings  are  accounted  for  through  the  introduction  of  the  Dirac  delta  or  impulse  function  as 
described  in  references  17  and  20.  The  cylinder  is  assumed  to  be  loaded  by  compressive 
and/or  pressure  loadings  (mechanical  loadii^),  and  is  subjected  to  thermal  stresses. 

In  considering  thermal  stresses  and  deformations,  the  usual  assumptions  of  engi¬ 
neering  thermal-stress  studies  are  adopted  in  that  coupling  between  heat  transfer  and 
deformation  is  ignored.  Thermal  effects  are  accounted  for  by  specifying  axisymmetric 
temperature  distributions  in  the  shell  and  stiffening  elements.  The  theory  is  sufficiently 
general  to  include  arbitrary  temperature  variations  through  the  thickness  of  the  shell  and 
depth  of  the  stiffening  elements.  The  effects  of  axisymmetric  temperature  variations 
along  the  length  of  the  shell  are  accounted  for  in  an  approximate  manner. 

The  theory  was  derived  by  specifying  strain  energy  expressions  corresponding  to 
nonlinear  Donnell-von  Karman  displacements  in  the  shell  and  stiffeners,  and  by  applying 
the  method  of  minimum  potential  energy  to  obtain  equilibrium  equations  and  associated 
boundary  conditions.  These  equations  were  separated  by  a  perturbation  procedure  into 
equations  governing  axisymmetric  behavior  prior  to  buckling  and  equations  governing 
behavior  at  buckling.  A  complete  development  of  the  theory  is  contained  in  appendix  B. 
Since  similar  developments  have  appeared  in  the  literature  (for  example,  refs.  17  and  21), 
only  pertinent  equations  are  presented. 

In  the  derivation  of  the  governing  equations,  the  ends  of  the  stiffened  cylinder  were 
considered  to  be  free  to  expand  longitudinally  (u  o).  With  this  assumption,  thermal 
buckling  is  a  consequence  of  circumferential  compressive  stresses  introduced  by  radial 
restraint  at  the  boundaries  (w  =  0)  and  restraints  resulting  from  diSerences  in  expansion 
between  the  stiffeners  and  the  shell. 
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Governing  Equations 


Prebuckling.-  The  equation  governing  the  prebuckling  behavior  of  a  multilayered 
stiffened  cylinder  subject  to  axisymmetric  thermal  and  mechanical  loadings  derived  from 
the  present  theory  (see  appendix  B  (eq.  (B21))  can  be  written  as 


2  (^11^12 

rI  Cii 


1 

R 


+ 


+  ^(Nx-Nts-Ntx) 


+  N-j'(x)  - 


P 


(1) 


The  equation  is  a  linear,  fourth-order,  ordinary  differential  equation  in  w^,  the  radial 
displacement  of  the  cylinder  prior  to  buckling.  The  coefficients  of  w^  involve  the 
structural  stiffnesses  of  the  composite  shell  wall  (including  stiffeners)  Dij,  Cy,  and 
Kij  which  are  associated  with  bending,  extension,  and  coupling  between  extension  and 
bending,  respectively.  Because  rings  are  treated  as  discrete  members,  the  circumfer¬ 
ential  extensional  stiffness  C22  is  a  function  of  x,  the  axial  coordinate.  Thus,  equa¬ 
tion  (1)  contains  involved  variable  coefficients  and  requires  numerical  techniques  for 
solution.  The  magnitude  of  the  applied  axial  compressive  load  appears  on  the  left- 
hand  side  of  equation  (1)  as  Nxw]^,  a  term  which  represents  the  contribution  of  nonlinear 
von  Karman  strain-displacement  terms  to  the  prebuckling  analysis. 

The  right-hand  side  of  equation  (1)  involves  loading  terms  resulting  from  tempera¬ 
ture  gradients  in  the  stiffened  shell  structure  (terms  with  subscript  T)  and  the  applied 
external  pressure  p.  The  thermal  loads  N-jy,  N>j>g,  and  N-j-j.  are  results  of 

temperature  gradients  through  the  thickness  of  the  shell  and  depth  of  the  stiffening  ele¬ 
ments,  whereas  Nt(x)  contains  terms  associated  with  exponential  variations  of  temper¬ 
ature  along  the  length  of  the  stiffened  shell.  The  effect  of  thermal  gradients  and  mechan¬ 
ical  loadings  is  both  to  stress  the  cylinder  and  change  its  shape  prior  to  buckling.  The 
total  effect  can  be  accounted  for  by  finding  w^  and  its  derivatives  by  solving  equa¬ 
tion  (1)  for  suitable  boundary  conditions  and  applied  loadings. 

Buckling.-  The  equations  governing  the  buckling  behavior  of  a  stiffened  shell 
derived  from  the  present  theory  (see  eqs.  (B22))  can  be  written  as 


^Xb,x  ^^yB,y  =  ® 

(2a) 

^^yB,x  ^yB,y  = 

(2b) 
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(2c) 


Equations  (2a)  and  (2b)  involve  the  stress  resultants  associated  with  extension  and 

Ny  \  and  shear  f^xygj  during  buckling.  Equation  (2c)  involves  the  moment  resultants 


and  contains  some  of  the  effects  of  deformations  and 


stresses  in  the  axisYuimetric  prebuckling  state  (terms  with  subscript  A).  In  form, 
equations  (2)  reduce  to  more  familiar  classical  buckling  equations  if  the  term  wV 


is  omitted  in  equation  (2c).  The  classical  buckling  theory  contains  the  effects  of  pre¬ 
buckling  loads  ^Nx,  Ny^>  ^^Ya)  prebuckling  deformations  /w^  and  w^j 


The  effects  of  thermal  gradients  appear  in  equations  (2)  in  two  ways.  First,  the 
buckling  stress  and  moment  resultants  are  functions  of  w^  so  that  terms  corresponding 
to  thermal  prebuckling  deformations  appear  in  equations  (2a),  (2b),  and  (2c).  Secondly,  a 
direct  thermal-stress  term  appears  as  Ny^Wg  in  equation  (2c).  This  term  repre¬ 
sents  the  effects  of  circumferential  stresses  introduced  by  differences  in  expansion  due  to 
thermal  gradients  or  variation  in  thermal  properties  in  the  stiffened  shell  wall.  In  many 
existing  theories  for  thermal  buckling,  the  direct  thermal-stress  term  is  retained, 
whereas  the  less  obvious  effects  of  thermal  deformations  are  omitted. 


Numerical  Solutions 

The  complexity  of  equations  (1)  and  (2)  requires  the  use  of  numerical  techniques. 
Solutions  to  the  equations  governing  prebuckling  as  well  as  buckling  behavior  were 
obtained  by  omitting  the  applied  torsional  loadings  l^xy^j  and  employing  finite  differences. 

Both  sets  of  equations  were  formulated  as  a  system  of  second-order  difference  equations 
which  are  documented  in  appendix  C.  The  prebuckling  equations  (eqs.  (B18)  and  (B21)) 
were  solved  by  matrix  algebra  by  exploiting  Gaussian  elimination.  Buckling  loads  were 
extracted  from  the  equations  governing  buckling  (eqs.  (C15)  and  (C16))  by  either  of  two 
methods:  determinant  plotting  or  modaT iteration.  In  the  determinant  plotting  solution, 
either  the  applied  axial  compressive  load  Nx,  the  pressure  p  or  a  base  temperature 
T^  can  be  used  as  a  buckling  parameter.  For  the  modal  iteration  solution,  T^  must 
be  the  buckling  parameter.  The  parameter  T^  is  related  to  temperature  changes 
between  the  cylinder  shell  and  stiffeners  and  the  boundaries  of  the  stiffened  cylinder. 
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The  two  numerical  solutions  have  been  incorporated  into  computer  programs 
entitled  BAMSOC  I  and  BAMSOC  II  (Buckling  Analysis  of  Multilayered  Stiffened 
Orthotropic  Cylinders).  The  programs  were  written  in  FORTRAN  IV  for  the  Control 
Data  6600  computer.  Input  variables  and  sample  problems  are  given  in  appendix  D.  The 
cards  and  computer  listings  for  these  programs  can  be  obtained  from  COSMIC,  University 
of  Georgia,  Athens,  Georgia  30601. 

APPLICATIONS  OF  ANALYSIS 

By  use  of  the  present  analysis,  the  thermal  buckling  behavior  of  four  types  of 
cylindrical  shells  was  investigated:  unstiffened,  ring  stiffened,  stringer  stiffened,  and 
ring  and  stringer  stiffened.  Unstiffened  and  ring-stiffened  shells  were  studied  in  order 
to  relate  results  from  the  present  theory  to  published  results  from  less  complex  theories. 
The  behavior  of  two  large -diameter  shells  —  one  stringer  stiffened,  the  other  ring  and 
stringer  stiffened  -  was  also  studied.  The  interaction  of  thermal  and  mechanical  loading 
in  these  two  structures  was  investigated  in  order  to  assess  the  sensitivity  of  contempo¬ 
rary  aerospace  designs  to  thermal  buckling. 

Unstiffened  Cylinders 

A  solution  for  buckling  of  unstiffened  cylinders  subjected  to  uniform  temperature 
has  been  presented  by  Hoff  in  reference  1.  In  the  reference,  calculations  of  the  thermal 
buckling  temperature  were  made  for  a  steel  cylinder  with  a  radius-thickness  ratio  of  300 
and  a  length-radius  ratio  of  about  0.3.  The  boundary  conditions  employed  corresponded 
to  those  of  simple  support  ^N^^  =  Vg  =  Wg  =  Mx^  =  0^ .  The  cylinder  was  found  to  buckle 

at  an  unusually  high  temperature. 

In  calculations  using  the  present  theory,  both  the  determinant  plotting  solution  and 
the  modal  iteration  solution  were  employed  in  an  attempt  to  find  a  buckling  temperature 
for  the  cylinder.  The  modal  iteration  solution  would  not  converge.  Further  investiga¬ 
tion  was  conducted  with  the  determinant  plotting  solution,  but  no  buckling  temperature 
could  be  found  within  a  practical  temperature  range. 

The  basic  difference  between  the  present  theory  and  that  of  reference  1  lies  in  the 
treatment  of  the  effects  of  deformations  and  stresses  in  the  shell  prior  to  buckling.  In 
the  theory  of  reference  1,  the  effects  of  thermal  stresses  introduced  by  restraints  offered 
by  the  boundaries  of  the  cylinder  are  accounted  for  by  introducing  circumferential  thermal 
stresses  into  Donnell's  buckling  equation.  In  the  Donnell  equation,  however,  the  effects  of 
rotations  and  changes  in  rotations  ^w^  and  w 

ignored.  The  present  theory  accounts  for  thermal  effects  in  a  more  consistent  fashion, 
and  as  previously  discussed,  contains  additional  terms  in  the  buckling  equations  which 


in  the  cylinder  prior  to  buckling  are 
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are  associated  with  shell  prebuckling  rotations.  In  order  to  make  computations  for  the 
steel  cylinder  on  a  more  comparable  basis,  these  additional  terms  were  suppressed  in 
buckling  calculations.  Thermal  buckling  predictions  were  then  in  excellent  agreement 
with  the  results  of  reference  1.  It  must  be  concluded  that  the  additional  terms  of  the 
present  theory  tend  to  alleviate  the  effects  of  the  boundary  restraints;  thus,  the  buckling 
temperature  of  this  cylinder  is  beyond  a  practical  range  of  interest. 


An  approximate  solution  for  buckling  of  clamped,  unstiffened  cylinders  subjected  to 
uniform  temperature  is  presented  by  Sunakawa  in  reference  6.  The  solution  is  based  on 
a  theory  which  includes  the  prebuckling  displacement  effect  just  discussed.  The  clamped 
boundary  conditions  for  buckling  implied  in  the  solution  are  =0,  =  0,  w^  =  0, 

B  15  -D 


and  =  0.  Thermal  buckling  predictions  for  the  steel  cylinder  previously  described 

were  performed  by  using  the  same  clamped  boundary  conditions  in  the  present  theory. 


Results  for  the  nondimensional  buckling  temperature  difference 
circumferential  waves  n  in  the  buckling  mode  are  as  follows: 


qT 


cr 


and  the  number  of 


Present 

theory 

Theory  of 
reference  6 

oTcr 

n 

0.01440 

46 

0.01255 

44 

It  can  be  seen  that  the  present  theory  yields  a  buckling  temperature  which  is  about  15  per¬ 
cent  higher  than  that  predicted  by  reference  6  with  similar  circumferential  mode  shapes. 
The  differences  in  buckling  predictions  are  attributed  to  lack  of  convergence  of  the  one- 
term  GalerkLn  solution  of  reference  6.  Convergence  of  the  present  solution  was  investi¬ 
gated  by  varying  the  number  of  finite-difference  stations  and  was  found  to  be  satisfactory 
when  100  finite-difference  stations  were  employed. 

Ring-Stiffened  Cylinders 

A  theoretical  model  for  buckling  of  ring- stiffened  cylinders  has  been  proposed  by 
Anderson  in  reference  10.  In  Anderson^s  theory,  the  ring-stiffened  cylinder  is  idealized 
by  considering  the  behavior  of  a  skin  bay  in  the  central  regions  of  an  infinitely  long  cyl¬ 
inder.  Thermal  stresses  caused  by  temperature  differences  between  the  skin  and  rings 
are  found  for  the  bay  by  assuming  that  no  rotation  occurs  at  the  rings.  The  circumfer¬ 
ential  stresses  are  then  introduced  into  DonnelPs  buckling  equation  for  an  unstiffened, 
isotropic  shell.  Boundary  conditions  (at  the  rings)  during  buckling  are  assumed  to  be 
those  corresponding  to  conventional  simply  supported  or  clamped  boundaries. 

In  contrast,  the  present  theory  treats  a  cylinder  of  finite  length.  Thermal  stresses 
and  deformations  prior  to  buckling  are  found  for  a  cylinder  supported  at  its  ends  in 
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conventional  fashion,  but  with  intermediate,  discrete  stiffening  rings.  Rotation  of  the 
rings  may  occur  during  thermal  loading.  Thermal  stresses  and  deformations  for  a  dis¬ 
cretely  stiffened  shell  are  introduced  consistently  into  the  buckling  equation.  Boundary 
conditions  during  buckling  are  enforced  only  at  the  ends  of  the  cylinder  so  that  the  con¬ 
tinuity  conditions  at  each  ring  are  natural  consequences  of  the  structural  stiffnesses  of 
the  ring  and  shell  at  the  line  of  attachment. 

Comparison  with  reference  10.-  In  order  to  compare  Anderson's  predictions  with 
the  present  results,  an  analytical  investigation  was  made  of  the  thermal  buckling  behavior 
of  ring- stiffened  cylinders  with  proportions  similar  to  the  test  cylinders  of  reference  11. 
These  cylinders  were  internally  stiffened  with  equally  spaced  Z-shaped  rings  having  the 
dimensions  shown  in  figure  2.  Geometric  and  mechanical  properties  employed  in  the 
present  calculations  are  indicated  in  table  I.  In  the  investigation,  the  rings  and  bound¬ 
aries  of  the  cylinder  were  considered  to  be  at  room  temperature,  whereas  the  cylinder 
skin  was  heated  to  a  constant  temperature.  The  nonuniform  temperature  in  the  stiffened 
cylinder  produces  circumferential  thermal  stresses  that  cause  buckling  of  the  cylinder 
skin  between  rings.  The  cylinder  was  assumed  to  be  simply  supported  at  its  ends  with 
the  boundary  conditions  Nxg  =  0,  Vg  =  0,  Wg  =  0,  and  Mxg  =  0. 


For  the  comparison,  the  length,  radius,  and  thickness  of  the  cylinder  were  held  con¬ 
stant  for  all  the  computations.  The  number  of  rings  in  the  cylinder  (and  hence  the  ring 

spacing  t)  was  varied  to  study  the  effects  of  the  bay  curvature  parameter 


on  thermal  buckling  characteristics.  The  results  of  reference  10  suggest  that  somewhat 
unusual  behavior  occurs  for  small  Z  in  ring-stiffened  cylinders.  The  overall  cylinder 
length  a  was  much  shorter  than  that  of  the  test  cylinders  of  reference  11.  The  particu¬ 
lar  value  selected  (see  fig.  2)  was  based  on  considerations  of  the  number  of  rings  and 
finite-difference  stations  required  to  investigate  small  values  of  Z. 


Results  of  the  theoretical  buckling  predictions  from  reference  10  and  the  present 
theory  are  compared  in  figure  3.  The  nondimensional  skin— buckling  temperature  oTcx* 
is  shown  as  a  function  of  Z.  Results  from  the  present  theory  are  shown  as  discrete 
points  corresponding  to  calculations  with  a  finite  number  of  rings.  The  solid  curve  has 
been  faired  through  the  points  to  indicate  trends.  The  dashed  curves  on  the  figure  are 
based  on  results  presented  by  Anderson  (ref.  10)  for  simply  supported  and  clamped  cylin¬ 
der  bays  with  rigid  rings.  A  correction  which  accounts  for  flexibility  of  the  rings  is  sug¬ 
gested  in  the  reference.  When  the  rings  are  considered  to  be  flexible,  circumferential 
thermal  stresses  near  the  rings  are  less  than  those  computed  for  rigid  rings;  therefore, 
the  thermal  buckling  temperatures  are  higher.  For  the  present  cylinders  the  correction 
amounts  to  an  increase  of  about  7  to  7^  percent  above  the  values  shown  in  figure  3. 


12 


In  figure  3 ,  it  can  be  seen  that  results  from  the  present  theory  are  in  reasonable 
agreement  with  Anderson's  results  for  clamped  cylinder  bays  when  Z  is  greater  than 
15.  The  numerical  differences  shown  for  this  range  of  Z  are  about  5  percent.  For 
smaller  Z,  a  larger  disparity  exists.  Studies  of  prebuckling  stress  distributions  based 
on  the  present  theory  suggest  that  the  idealizations  used  by  Anderson  to  obtain  thermal 
stresses  are  reasonably  accurate.  It  is  believed,  therefore,  that  the  analytical  differ¬ 
ences  shown  in  figure  3  stem  mainly  from  the  additional  thermal  deformation  terms  in 
the  present  buckling  theory,  and  the  effective  ring  boundary  conditions  during  buckling. 

The  present  results  suggest  that  when  Z  is  greater  than  l5,  the  cylinder  rings  offer 
nearly  clamped  support  to  the  skin  bays  during  buckling.  Furthermore,  both  the  present 
theory  and  that  of  reference  10  suggest  that  the  thermal  buckling  temperature  becomes 
independent  of  ring  spacing  for  this  range  of  Z.  The  theoretical  differences  and  unusual 
behavior  for  cylinders  with  smaller  Z  is  probably  associated  with  the  interaction  of 
closely  spaced  rings. 

Thermal-mechanical  interactions.  -  A  further  topic  of  interest  in  the  thermal  buck¬ 
ling  behavior  of  ring- stiffened  cylinders  is  the  interaction  between  thermal  and  mechan¬ 
ical  loading.  The  behavior  pf  the  simply  supported  cylinders  just  discussed  was  investi¬ 
gated  under  combinations  of  axial  compression  and  thermal  loads  induced  by  heating  the 
cylinder  skin.  Results  of  interaction  studies  are  presented  in  figure  4,  where  the  nondi- 
mensional  mechanical  buckling  load  Nx^^^Et  has  been  plotted  against  the  thermal  buck¬ 
ling  load  oTcr*  Buckling  interaction  curves  are  presented  for  cylinders  of  fixed  overall 
geometry  (see  fig.  2)  with  three,  five,  and  seven  rings.  The  buckling  loads  shown  for 
axial  compression  loading  alone  =  oj  were  compared  with  published  results  for  iso¬ 

tropic  cylinders  having  a  similar  value  of  Z  (ref.  22).  The  comparison  suggested  that 
the  rings  of  the  present  cylinders  provide  almost  clamped  supports  for  the  skin  bays 
when  loaded  in  axial  compression.  A  similar  result  has  already  been  noted  for  thermal 
loading  when  Z  is  large.  The  interaction  curves  shown  are  very  regular  in  shape  with 
the  exception  of  the  case  of  three  rings.  Results  of  the  interaction  curves  were  compared 
with  interactions  suggested  by  reference  10  for  clamped  cylinders  but  no  correlation  was 
evident  for  the  range  of  curvature  parameter  considered  (Z  =  8  to  30).  In  this  range  of 
Z,  the  results  of  reference  10  suggest  that  substantially  larger  thermal  loads  can  be  car¬ 
ried  when  axial  compressive  loads  are  present,  whereas  results  from  the  present  theory 
do  not  support  this  conclusion. 


Stringer-Stiffened  Cylinder 

In  order  to  investigate  the  thermal  buckling  behavior  of  longitudinally  (stringer-) 
stiffened  shells,  the  behavior  of  a  large-diameter  aluminum  shell  representative  of  pro¬ 
portions  which  might  be  considered  for  an  interstage  structure  in  a  large  launch  vehicle 
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was  investigated.  The  structure  considered  is  an  integrally  stiffened  shell  which  has 
been  used  in  theoretical  and  experimental  studies  of  stiffener  and  loading  eccentricity 
effects.  (See  refs.  23  and  24.)  The  geometry  and  mechanical  properties  adopted  for  the 
present  investigation  are  given  in  figure  5  and  table  1(b).  Thermal  buckling  problems  in 
interstage  structures  might  conceivably  arise  from  thermal  gradients  induced  by  rapid 
aerodynamic  heating  of  the  shell  outer  skin;  as  a  result,  thermal  lags  between  stiffener 
and  shell  temperatures  occur.  For  the  present  analysis,  the  skin  was  assumed  to  be 
heated  whereas  the  stringers  remained  unheated.  Bending  and  thrust  loads  are  also 
present  during  flight;  thus,  it  seemed  desirable  to  investigate  the  interaction  between 
axial  compressive  and  thermal  loads.  For  purposes  of  the  investigation,  the  cylinder 
was  assumed  to  be  clamped  at  its  ends  ^Nxg  =  Vg  =  Wg  =  Wg^  =  oj.  Local  buckling  of 

the  skin  was  not  considered. 

Prebu^kling  stress  distribution.-  The  stress  distributions  induced  by  combinations 
of  axial  compression  and  temperature  were  investigated  by  use  of  the  results  of  the  solu¬ 
tion  of  the  equations  governing  deformations  and  stresses  prior  to  buckling.  (See  eq.  (1).) 
Stress  distributions  in  the  skin  and  stringers  of  the  interstage  for  three  combinations  of 
applied  compressive  load  Njj  and  change  in  skin  temperature  Tgjj  are  presented  in 
figure  6  as  functions  of  the  distance  x/a  measured  from  the  supported  end  of  the  shell. 
In  figure  6(a)  the  axial  stress  at  the  skin  middle  surface  is  shown  to  increase  in  propor¬ 
tion  to  the  applied  axial  load.  The  axial  stress  asymptotically  approaches  the  mem¬ 
brane  state  of  stress  =  Nx/Et)  in  the  central  areas  of  the  shell.  The  higher  values  of 
stress  near  the  ends  of  the  shell  (x/a  =  0)  are  due  to  bending  stresses  developed  from  end 
constraints.  In  figure  6(b),  the  circumferential  stresses  at  the  skin  middle  surface  shown 
increase  roughly  in  proportion  to  skin-temperature  changes.  The  circumferential 
stresses  are  small  in  the  central  regions  of  the  cylinder  and  the  stresses  shown  near  the 
ends  are  largely  due  to  restraints  in  expansion  between  the  cold  boundaries  and  cold 
stringers  and  the  heated  shell.  A  significant  difference  between  circumferential  stress 
variation  in  the  present  stiffened  shell  and  results  for  unstiffened  shells  of  similar  geom¬ 
etry  is  that  stresses  decay  much  more  rapidly  in  unstiffened  shells.  (See  ref.  13.) 

Stresses  developed  in  the  longitudinal  stiffening  for  various  combinations  of  axial 
compression  and  thermal  load  are  shown  in  figure  6(c).  The  stress  distribution  in  the 
stringer  is  plotted  against  the  normal  coordinate  z  nondimensionalized  by  the  stringer 
height  h  for  various  longitudinal  stations  x/a.  The  origin  of  z  was  taken  to  be  the 
surface  of  the  stiffener  in  contact  with  the  shell.  The  stress  distributions  shown  in  fig¬ 
ure  6(c)  indicate  that  bending  stresses  in  the  shell  increase  appreciably  as  the  tempera¬ 
ture  difference  between  skin  and  stringers  increases.  For  the  case  with  no  applied  axial 
compression,  the  calculated  stress  in  the  outstanding  flange  of  the  stiffener  actually 
exceeds  the  ultimate  tensile  strength  of  aluminum. 
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The  results  presented  for  the  stress  distribution  in  the  interstage  structure  suggest 
that  complex  states  of  stress  exist  which  depend  on  the  relative  magnitude  of  the  applied 
axial  compression  and  the  differences  in  expansion  between  the  boundaries  and  stringers, 
and  the  heated  shell.  Studies  of  the  detailed  stress  distribution  in  skin  and  stiffeners 
indicate  that  intolerably  large  bending  stresses  exist  in  the  stringers  near  the  ends  when 
only  thermal  loads  are  present.  The  results  also  indicate  that  thermal  stresses  are 
developed  from  boundary  effects  which  appear  to  persist  over  relatively  large  sections 
of  the  longitudinally  stiffened  shell. 

Buckling  of  stringer -stiffened  cylinder.-  The  buckUhg  behavior  of  the  idealized 
interstage  structure  was  investigated  under  combinations  of  axial  compressive  and  ther¬ 
mal  loads.  Results  of  the  investigation  are  presented  in  figure  7  as  an  interaction  curve 
between  the  normalized  axial  compressive  load  at  buckling  change  in 

skin  temperature  at  buckling  oTcj..  The  three  combinations  of  axial  load  and  tempera¬ 
ture  change  discussed  in  the  section  entitled  "Prebuckling  stress  distribution"  are  shown 
as  cross  marks  in  figure  7.  The  buckling  load  shown  for  the  clamped  shell  under  axial 
compression  alone  »  o)  was  in  good  agreement  with  published  results,  the  present 

calculation  being  about  2  percent  higher  than  the  result  presented  in  reference  24. 

In  figure  7  it  can  be  seen  that  the  interaction  curve  has  an  unusual  shape  when  small 
thermal  loads  are  present.  For  temperature  changes  of  about  50°  F  (28  K)  in  the  skin, 
that  is,  oTcr  =  0.0007,  the  shell  buckling  load  actually  increases  to  about  1.4  times  the 
compressive  buckling  load  of  an  unheated  shell.  It  appears  that  small  thermal  loads  can 
introduce  favorable  deformations  prior  to  buckling  which  increase  the  shell’s  effective 
stiffness  to  resist  axial  compression.  As  can  be  seen  from  the  interaction  curve,  how¬ 
ever,  when  thermal  loads  became  large,  the  longitudinally  stiffened  shell  had  little 
resistance  to  the  circumferential  thermal  stresses  so  that  buckling  can  occur  by  thermal 

=  o| .  The  temperature  change  corresponding  to  this  case  is  about 

130°  F  (72  K),  and  is  well  within  a  practical  temperature  range  for  aluminum.  In  inter¬ 
stage  structures,  when  thermal  environments  are  present,  simple  insulators  such  as  cork 
are  used  to  shield  the  outer  skin  of  the  structure.  However,  in  many  instances,  the  insu¬ 
lation  is  not  adequate  to  eliminate  all  thermal  gradients.  The  present  results  suggest 
that  small  gradients  may  be  tolerated,  but  if  large  thermal  gradients  are  present,  the 
structure  should  be  carefully  considered  with  regard  to  the  possibility  of  thermal  buckling. 

The  normal  component  of  the  buckling  mode  shape  for  the  interstage  structure  under 
thermal  loading  alone  is  shown  in  figure  8,  where  the  normalized  variation  of  the  radial 
buckling  displacements  ^/^raax  plotted  against  the  length  parameter  x/a  where 
a  is  the  length  of  the  shell.  For  comparison,  prebuckling  deformations  at  the  onset  of 
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buckling  are  shown  in  the  inset  in  figure  8.  The  buckling  wave  pattern  is  characterized 
by  a  nearly  sinusoidal  wave  along  the  length  of  the  shell  and  many  sinusoidal  waves 
around  the  circumference  (n  =  65).  Because  the  number  of  circumferential  waves  is 
large,  the  nodes  in  the  buckle  pattern  in  the  circumferential  direction  are  separated  by 
distances  which  are  about  times  the  stringer  spacing.  Since  the  stiffness  properties 
of  longitudinal  stiffening  are  averaged  in  the  present  theory,  the  calculated  circumferen¬ 
tial  mode  shape  is  one  which  is  near  the  limits  of  applicability  of  the  theory. 

Ring-  and  Stringer -Stiffened  Cylinder 

In  order  to  investigate  the  thermal  buckling  behavior  of  ring-  and  stringer-stiffened 
shells,  the  behavior  of  a  stiffened  cylinder  with  geometry  similar  to  that  which  might  be 
considered  in  a  preliminary  design  of  a  supersonic  transport  was  investigated.  A  rather 
long,  uniform,  titanium  fuselage  section,  bounded  forward  by  a  heavy  bulkhead  and  aft  by 
the  wing-joint  structure,  was  selected.  The  geometry  and  mechanical  properties  adopted 
for  the  present  calculations  are  presented  in  figure  9  and  table  I.  Studies  of  design 
loads  and  predicted  thermal  environment  for  this  type  of  structure  suggested  the  possi¬ 
bility  of  thermal  buckling  as  a  consequence  of  a  rapid  climb  to  altitude.  In  this  maneuver 
there  is  a  thermal  lag  between  the  skin  and  stiffening  elements  so  that  for  a  short  period 
of  time  the  skin  is  at  elevated  temperatures  while  both  rings  and  stringers  are  unheated. 
Bending  and  pressurization  loads  on  the  fuselage  are  also  present  during  the  maneuver; 
therefore,  it  seemed  desirable  to  investigate  the  interaction  of  axial  compressive  loads 
and  thermal  loads  for  various  internal  pressures.  For  purposes  of  the  investigation,  the 
cylinder  was  assumed  to  be  simply  supported  at  its  ends  ^Nxg  =  Vg  =  '^■q  =  ^ 

the  axial  compression  was  introduced  by  means  of  a  line  load  applied  at  the  skin  middle 
surface. 

Prebuclding  stress  distribution.-  The  stress  distribution  induced  by  combinations 
of  axial  compression  and  temperature  was  investigated  by  using  results  based  on  the 
solution  of  equation  (1).  Stress  distributions  in  the  skin,  stringers,  and  rings  of  the  fuse¬ 
lage  under  an  applied  compressive  load  of  1650  Ibf/in.  (0.289  MN/m)  and  temperature 
difference  of  360°  F  (200  K)  between  skin  and  stiffeners  were  obtained.  Results  are  pre¬ 
sented  in  figure  10.  The  diffusion  of  load  into  each  element  of  the  structure  is  illustrated 
by  showing  results  for  various  stations  near  one  of  the  supported  ends  of  the  cylinder. 

In  figure  10(a)  the  stress  distribution  in  the  skin  of  the  cylinder  is  shown  as  a  function  of 
length.  At  the  end  of  the  cylinder  (x/a  =  0),  both  axial  and  circumferential  stresses  are 
very  high,  but  are  still  within  the  elastic  range  of  titanium.  The  stresses  are  not  nearly 
as  large  as  those  for  the  clamped  interstage  structure  already  discussed.  Sharp  peaks  in 
stress  occur  at  each  of  the  19  rings  in  the  shell.  The  circumferential  stress  shown  in 
bays  that  are  distant  from  the  boundary  is  a  result  of  the  differences  in  thermal  expansion 
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between  the  rings  and  the  shell.  The  curves  shown  in  figure  10(a)  indicate  that  the  state 
of  stress  does  not  vary  from  bay  to  bay  beyond  the  third  bay  (x/a  *  0. 15). 

Variations  in  the  stress  distribution  in  the  stiffening  elements  are  shown  in  fig¬ 
ure  10(b).  The  assumptions  of  the  present  theory  are  such  that  the  ring  is  in  a  state  of 
uniform,  uniaxial  stress  (in  this  case,  tension).  The  ring  stress  variations  shown  in  fig¬ 
ure  10(b)  indicate  that  the  interior  ring  stresses  are  essentially  constant  in  rings  that  are 
not  adjacent  to  the  boundary.  The  assumptions  of  the  present  theory  are  such  that  the 
stringers  (longitudinal  stiffening)  are  also  in  a  state  of  uniaxial  stress,  although  bending 
stresses  through  the  depth  of  the  stringers  may  exist.  The  magnitudes  of  these  bending 
stresses  are  given  in  figure  10(b)  by  plotting  stresses  against  the  depth  parameter  z/h 
for  various  longitudinal  stations.  The  bending  stresses  in  the  interior  of  the  shell  are 
mainly  due  to  thermal  stresses  and  restraints  to  expansion  offered  by  the  rings.  The 
bending  stresses  near  the  ends  of  the  cylinder  also  include  the  effects  of  introducing  the 
applied  axial  compression  load  by  means  of  a  line  load  on  the  skin  midplane.  The  edge 
moment  introduced  by  the  eccentricity  of  the  line  load  with  respect  to  the  stringer-shell 
centroid  has  been  noted  in  previous  studies  (ref.  23)  and  large  effects  on  buckling  predic¬ 
tions  have  been  found.  In  the  present  case,  the  stress  distributions  shown  in  figure  10(b) 
suggest  that  this  edge  moment  decays  rapidly  from  the  cylinder  ends  so  that  the  state  of 
stress  in  the  stringers  does  not  vary  from  bay  to  bay  beyond  the  second  ring. 

The  results  presented  for  the  stress  distribution  in  this  example  indicate  the  com¬ 
plexity  of  the  state  of  stress  existing  in  a  ring-  and  stringer-stiffened  shell  under  com¬ 
bined  mechanical  and  thermal  loadings.  Studies  of  the  detailed  stress  distributions  in 
skin  and  stiffening  elements  suggest  that  no  large  stresses  exist  that  are  not  within  the 
capabilities  of  the  structural  material.  The  effects  of  the  boundaries  of  the  stiffened 
shell  on  stress  distribution  appear  to  be  negligible  beyond  the  second  bay  of  the  cylinder. 

Buckling  of  ring-  and  stringer-stiffened  cylinder.-  The  buckling  behavior  of  the 
idealized  supersonic  transport  design  was  investigated  under  combinations  of  axial  com¬ 
pressive  and  thermal  loads.  In  addition,  the  effects  of  including  cabin  pressurization  as 
an  applied  load  were  also  included.  The  results  of  the  investigation  are  summarized  in 
figure  11  by  means  of  interaction  curves  between  axial  compressive  load  at  buckling 
Nx(.j.|Ef  and  relative  skin  temperature  at  buckling  oT^r-  The  solid  curves  shown  cor¬ 
respond  to  various  values  of  internal  pressure  p  normalized  by  sea-level  atmospheric 
pressure  p^.  The  maximum  value  of  p/p^  shown  is  about  that  proposed  for  cruising 
altitude. 

Results  shown  in  figure  11  for  p/Pq  =  0  indicate  that  the  Interaction  between  axial 
compressive  and  thermal  loads  is  substantial.  In  fact,  in  the  absence  of  axial  compres¬ 
sion,  thermal  loads  are  large  enough  to  cause  buckling  in  a  practical  temperature  range. 
The  value  of  oTgr  shown  for  this  case  corresponds  to  a  temperature  difference  of  about 
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525°  F  (292  K)  between  the  skin  and  stiffening  elements.  Results  shown  in  figure  11  for 
cases  where  internal  pressurization  exists  ^p/Pq  o)  indicate  that  even  small  values  of 
internal  pressure  are  very  beneficial  in  alleviating  the  problem  of  buckling  under  thermal 
loading  alone.  At  cabin  pressure  at  cruise  =  -0.82),  the  interaction  curve  is  nearly 

horizontal  over  a  practical  temperature  range  and  the  buckling  load  is  roughly  1-  times 

Li 

that  of  an  unpressurized  fuselage.  In  some  aircraft  design  criteria  for  buckling,  how¬ 
ever,  the  beneficial  effects  of  pressure  are  not  included  in  buckling  calculations;  there¬ 
fore,  the  fuselage  structure  is  required  to  be  stable  even  if  all  cabin  pressure  is  acciden¬ 
tally  lost.  Under  these  circumstances,  elimination  of  buckling  under  thermal  and  com¬ 
pressive  loads  or  thermal  loads  alone  may  be  a  serious  design  consideration. 

An  additional  effect  shown  in  figure  11  by  the  dashed  curve  is  that  of  introducing  the 
axial  compressive  load  at  the  neutral  surface  of  the  sheet-stringer  combination 
(e  =  0.13  in.  (0.33  cm))  rather  than  at  the  skin  midplane  (e  =  0).  The  dashed  curve  sug¬ 
gests  substantial  improvement  in  buckling  strength  under  axial  compression  alone  due  to 
the  edge  moment  effect  already  mentioned.  However,  the  interaction  curve  is  unaffected 
when  large  thermal  loadings  are  present. 

The  normal  component  of  the  buckling  mode  shape  for  the  stiffened  shell  under 
thermal  loading  alone  with  no  internal  pressure  is  presented  in  figure  12,  where  the  nor¬ 
malized  variation  of  the  radial  buckling  displacements  ^^Wnlax  plotted  against 
length  x/a  for  one-half  of  the  cylinder.  For  comparison,  prebuckling  deformations  at 
the  onset  of  buckling  are  shown  in  the  inset  in  figure  12.  The  wave  shape  shown  in  fig¬ 
ure  12  is  typical  of  a  panel  instability  buckling  pattern  in  that  nodes  occur  at  each  of  the 
rings.  It  can  be  seen  that  the  amplitude  of  W  decays  in  the  interior  bays  so  that  in  the 
central  bays  of  the  cylinder  almost  negligible  deformation  is  present.  The  variation  in 
amplitude  from  bay  to  bay  suggests  that  buckling  is  a  consequence  of  restraints  at  or 
near  the  boundaries  of  the  cylinder.  The  presence  of  longitudinal  stiffening  apparently 
makes  thermal  edge  effects  persist  over  a  large  length  as  evidenced  by  the  large  ampli¬ 
tude  of  the  buckle  wave  in  the  fifth  bay  of  the  cylinder  (0.20  ax/a  «  0.25).  As  was  the 
case  with  the  interstage  structure,  a  large  number  of  circumferential  waves  (n  =  57) 
occurred  in  the  buckling  mode  when  thermal  loading  predominated.  In  this  case,  however, 
at  least  three  stringers  are  included  in  each  circumferential  half-wave.  Thus,  the  buck¬ 
ling  mode  found  appears  to  be  a  legitimate  panel  instability  mode  and  is  consistent  with 
theory  employing  the  average  stiffness  properties  of  longitudinal  stiffening. 

The  present  study  of  an  idealized  supersonic  transport  design  indicates  that  a  poten¬ 
tial  design  problem  exists  in  that  temperature  gradients  which  might  be  induced  by  a 
rapid  climb  to  altitude  were  sufficient  to  buckle  the  shell  in  a  practical  temperature  range. 
The  presence  of  internal  cabin  pressure  is  sufficient  to  alleviate  the  thermal  loads.  If 
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internal  cabin  pressure  is  not  available  as  a  stabilizing  force,  additional  reinforcement 
must  be  judiciously  applied  to  the  structure. 

CONCLUDING  REMARKS 

An  extension  of  the  theory  of  NASA  Technical  Note  D-4283  for  buckling  of  stiffened 
cylindrical  shells  under  combinations  of  mechanical  and  thermal  loadings  has  been  pre¬ 
sented.  It  includes  the  effects  of  eccentricity  in  multilayered  reinforcement,  deforma¬ 
tions  prior  to  buckling,  and  discrete  rings.  The  theory  is  sufficiently  general  to  account 
for  the  averaged  stiffnesses  of  stringers  and  multiple  orthotropic  shell  layers  as  well  as 
arbitrary  temperature  distributions  through  the  thickness  of  the  shell  and  depth  of  the 
stiffeners.  Buckling  solutions  have  been  obtained  by  employing  finite  differences.  Two 
computer  programs  from  which  buckling  loads  can  be  found  (either  by  determinant  plotting 
or  modal  iteration)  have  been  developed. 

Theoretical  results  for  thermal  buckling  have  been  compared  with  existing  solutions 
for  unstiffened  and  ring-stiffened  cylindrical  shells.  Results  from  the  present  theory  for 
a  simply  supported  unstiffened  shell  indicated  that  thermal  buckling  would  not  occur  in  a 
practical  temperature  range.  Results  for  a  clamped  unstiffened  shell  were  in  agreement 
with  published  theoretical  results.  Buckling  predictions  for  a  ring- stiffened  shell  with 
various  ring  spacings  were  in  agreement  with  existing  theoretical  results  when  the  cylin¬ 
der  bay  curvature  parameter  is  large.  In  cases  where  large  disparities  between  pub¬ 
lished  theoretical  results  and  those  presented  herein  exist,  differences  are  attributed  to 
the  effects  of  prebuckling  deformations  which  have  not  been  accounted  for  in  previous 
theories. 

The  effects  of  interaction  of  thermal  loading  and  axial  compression  on  the  buckling 
behavior  of  two  types  of  contemporary  stiffened  shell  structures  have  been  investigated. 
Thermal  loads  were  introduced  by  considering  the  skin  of  the  structure  to  be  heated, 
while  the  stiffening  elements  and  shell  boundaries  remain  unheated.  The  behavior  of  an 
aluminum,  clamped,  large-diameter  longitudinally  stiffened  cylinder  representative  of  a 
launch  vehicle  interstage  structure  was  investigated.  Stress  distribution  studies  indi¬ 
cated  that  thermal  stresses  induced  by  moderate  skin  temperatures  cause  large  bending 
stresses  in  the  stiffeners  near  the  ends.  Buckling  calculations  suggest  that  the  inter¬ 
stage  structure  may  buckle  under  thermal  loading  alone  at  a  realistic  skin  temperature. 

The  behavior  of  a  titanium  simply  supported,  ring-  and  stringer-stiffened  cylinder 
representative  of  a  preliminary  design  for  a  supersonic  transport  fuselage  was  also 
investigated.  Stress  distribution  studies  indicated  that  no  large  thermal  stresses  beyond 
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the  capabilities  of  the  shell  material  were  present.  Buckling  calculations  suggest  that 
the  fuselage  also  can  buckle  under  thermal  loading  alone  at  a  practical  skin  temperature. 

Near  the  completion  of  this  work,  preliminary  calculations  for  the  ring-stiffened 
cylinders  by  another  investigator  suggested  that  thermal-buckling  behavior  could  be 
influenced  by  the  type  of  ring  theory  employed  in  the  analysis.  For  the  ring-stiffened 
cylinders,  out-of-plane  bending  stiffness  appeared  to  be  an  important  parameter.  To 
account  for  this  effect,  a  ring  theory  which  includes  out-of-plane  bending  behavior  was 
developed  and  is  presented  in  appendix  E.  This  ring  theory  was  incorporated  into  one  of 
the  computer  programs;  however,  no  numerical  results  were  available  for  this  paper. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Hampton,  Va.,  April  28,  1970. 
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APPENDIX  A 


CONVERSION  OF  U.S.  CUSTOMARY  UNITS  TO  SI  UNITS 

The  International  System  of  Units  (SI)  was  adopted  by  the  Eleventh  General  Con¬ 
ference  on  Weights  and  Measures  in  1960.  (See  ref.  18.)  Conversion  factors  for  the 
units  used  in  this  report  are  given  in  the  following  table: 


Physical  quantity 

U.S.  Customary 
Unit 

Conversion  factor 
{*) 

SI  Unit 
(**) 

Length . 

in. 

0.0254 

meters  (m) 

Stress  modulus . 

ksi 

6.895  X  106 

newtons/meter  2  (N/m2) 

Stress  resultant . 

Ibf/in. 

175.1 

newtons / meter  (N /m) 

Temperature  change  .... 

op 

5/9 

Kelvin  (K) 

^Multiply  value  given  in  U.S.  Customary  Unit  by  conversion  factor  to  obtain 


equivalent  value  in  SI  Units. 

**Prefixes  to  indicate  multiple  of  units  are  as  follows: 


Prefix 

Multiple 

giga  (G) 

109 

mega  (M) 

106 

kilo  (k) 

103 

deci  (d) 

10-1 

centi  (c) 

10-2 

milli  (m) 

10-3 
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APPENDIX  B 


DEVELOPMENT  OF  GOVERNING  EQUATIONS 


Stress-Strain  Relations 

In  the  ith  layer  of  the  shell,  the  stress-strain  relations  can  be  written  as 


1  -  u  u 


^y=— 

1  - 


7^  =G^  y* 

xy  xy'^xy  J 

where  cx^  and  oy  are  coefficients  of  linear  expansion  associated  with  thermal  elonga¬ 
tion  in  the  axial  and  circumferential  directions,  respectively,  T  is  a  change  in  temper¬ 
ature  relative  to  a  thermal-stress -free  state  and  the  remaining  quantities  with  the  super¬ 
script  i  are  orthotropic  elastic  constants  associated  with  extension.  The  strains 
appearing  in  equation  (Bl)  are  related  to  the  shell  displacements  u,  v,  and  w  by 
Donnell-von  Karman  relations  having  the  form 


e  =  ey  +  Z/Cy 


^xy  ^xy 


=  “.y  *  ''.X  +  »,x»,y 


%  ^,XX 
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APPENDIX  B  -  Continued 


'^y  =  -^,yy 


%y  - 


In  the  stringers 


whereas  in  the  rings 


%  = 


(B3) 


(B4) 


To  account  for  axisymmetric  temperature  variation  along  the  length  of  the  shell,  the 
relative  temperature  of  the  shell  is  expressed  as 


T  =  Tsk  1  +  ygj^e 


gx/al 


(B5) 


where  a  is  the  length  of  the  shell,  and  and  g  are  constants.  Similarly,  in  the 

stringers. 


T  =  Tg^  +  y  j 


(B6) 


and  in  the  rings 


T  =  Tr  1  +  Yr^ 


(B7) 


where  Tgj^,  Tg,  and  Tr  are  functions  of  z.  By  specifying  temperature  distributions 
with  equations  (B5)  to  (B7),  the  first-order  effects  of  longitudinal  variation  in  thermal 
distributions  can  be  obtained.  The  present  theory,  however,  does  not  consider  the  effects 
of  variation  of  mechanical  properties  with  the  assumed  longitudinal  temperature 
distribution. 


Potential  Energy 

The  integration  of  the  familiar  strain  energy  expression 


Hs 


1 

2 


roll 

Volume  L 


Ex 


+ 


.t,t 


+ 


_3[ — 

1  -  MxMyVy^ 


dV 
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APPENDIX  B  -  Continued 


over  the  depth  of  the  composite  shell  wall  (shell  plus  stiffeners)  yields  the  following 


expression: 


hs  =  I  j'y  +  NyCy  +  Nxyrxy  + 


Wx  +  Ny€y  +  Nxyrxy  +  ^x^x  +  ^y/Cy  +  (M^y  "  My^)^  "  (^x  +  ^gj^x  +  (^Tx  +  ^Ts^x 

Surface  v. 

^Nxy  +  Nxr^ey  +  ^Mxy  +  Mxrj  dx  dy  (gg) 


where 


“  ^11^  ^12^y  ■*■  ^11%  ^12^y  "  \^Tx  +  ^Ts 
Ny  =  Ci2ex  +  ^22^y  +  ^12«x  +  ^22i^y  '  f^Ty  + 


Nxy  =  Cee^'xy  +  ^66%y 

^X  =  Kll^  +  Ki2€y  +  DllKx  +  Dl2Ky  "  (M-Jx  +  M^g) 

My  =  Ki2ex  +  *^22ey  +  Dl2«x  +  D22«y  "  (^Ty  +  ^Tr] 


^xy  "^66^xy  ^GB'^xy  "*■  2d 

N 

Myx  =  ^66>'xy  ■'■  ^G&'^xy  ^  ^rJr%y 

j=l 

where  6(x  -  jZ)  is  the  Dirac  delta  function  associated  with  the  N  rings  on  the  cylinder. 

The  quantities  Cj^j,  D^j,  and  K^j  appearing  in  equations  (B9)  are  structural  stiff¬ 
nesses  associated  with  extension,  bending,  and  extension-bending  coupling,  respectively, 
and  are  given  by  the  following  integrals: 


1  .  ^S’^s 

Cll=  \-; - 6Z+-— 

I  1  “  MxMy 


dZ 


f^xMy 


Coo  =  - - - dz  +  )  6(x  -  jZ) 

J  1  jti 


Er  A-r 


(BlOa) 


(BlOb) 


(BlOc) 


*^66  "  J  ®xy  ^ 


(BlOd) 


APPENDIX  B  -  Continued 


Ex  EoAcZ, 

Kii  =  \ - ^ - z  dz  + 

“  ''-fXxiUy  d 

=  I  - z  dz 

1  -  ix^iXy 


N 


K22  =  — - - z  dz  +  y  6(x  -  jZ)  EfArZ. 

j  1 

*^66  =  J  ^xy^ 

E 


Dii  = 


X  _  ^2(jz  ■ 


1  -  fi^My 


d  d 


Di2  =  I  — ^ —  z^dz 
^  “  l^x^y 


E, 


N 


D22  =  1 - - - z^dz  +  y  6(x  -  jZ)|ErIr  +  EyArZ 

1  -  fXx/Xy  \ 


D66=y 


GjjyZ^dz 


(BlOe) 

(BlOf) 

(BlOg) 

(BlOh) 

(BlOi) 

(BlOj) 

(BlOk) 

(BIOI) 


where  the  integration  is  to  be  performed  over  the  multiple  layers  of  the  shell  wall.  A 
convenient  numerical  scheme  for  the  integration  for  a  specified  reference  surface  can  be 
found  in  equations  (83)  and  (85)  of  reference  19  (p.  34).  The  properties  associated  with 
the  stiffeners  (location  of  the  centroid  of  the  stiffener  z,  the  moment  of  inertia  I,  the 
cross-sectional  area  A,  and  the  torsional  constant  J)  are  also  defined  (when  appropri¬ 
ate)  with  respect  to  an  arbitrary  reference  surface.  In  specifying  the  structural  stiff¬ 
nesses,  the  dependence  of  the  elastic  constants  upon  temperature  can  be  accounted  for  in 
the  shell  by  analytically  dividing  the  multilayered  shell  into  thin  layers  and  performing 
the  required  integration  of  equations  (BIO)  numerically.  In  the  stiffeners,  the  following 
sample  definitions  can  be  employed  to  perform  a  comparable  integration  through  the  depth 
of  the  stiffeners: 
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APPENDIX  B  -  Continued 


£2^.  if  EsWdAs 

d  d  JAs 

Es(z)zclAs 

d  d  Jas 


The  thermal  moment  M.j.  and  thermal  force  terms  N.p  appearing  in  equa¬ 
tions  (B9)  are 

Ntx  = 

Nxy  =  CTy/l  + 

/  \  > 

%s  =  ^Ts(|^  +  ys^^J 

N  V 

%r  =  2  ^Trfl  + 

i=l  ^  J 


Mtx  - 

M.py  = 

%s  = 

Mfjip  = 


(B12) 
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For  Cfjijj,  C-py,  Krpjj,  and  Kj-y,  the  integration  is  to  be  performed  over  the 
multiple  layers  of  the  shell  wall  and  can  be  accomplished  by  the  method  suggested  for 
equations  (BIO). 

As  in  reference  17  the  potential  energy  of  mechanical  loading  for  the  case  of  an 

^  A 

applied  compression  load  Nx,  an  applied  shear  load  Nxy,  and  an  external  pressure  p 
is  given  by 

p2irR  pa  p27rR  ^  ,  v  a  p2TrR  ..  a 

”l  -  Jd  P*  *  +  Jo  0  ^  Jo  0  ® 

where  e  is  the  distance  from  the  arbitrary  reference  surface  to  the  line  on  which  the 
applied  mechanical  load  resultant  Nx  acts. 
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APPENDIX  B  -  Continued 


Equilibrium  Equations  and  Boundary  Conditions 


By  the  principle  of  minimum  potential  energy,  the  vanishing  of  the  first  variation  of 
the  total  potential  energy  (5(^3  +  IIl)  =  o)  yields  the  equilibrium  equations  and  appropriate 
boundary  conditions.  The  first  variation  with  respect  to  the  displacement  variables  u, 

V,  and  w  yields 


Nx,x  +  Nxy,y  -  ^ 


Nxy,x  +  Ny,y  -  0 


)  (B14) 


•M 


,xx  +  (^xy  -  Myx),xy  ‘  ^y,yy  "  NxW^xx  -  %^,yy  "  ^NxyW^xy  +  P  -  0 


with  the  following  boundary  conditions  to  be  specified  at  each  end  of  the  cylinder: 

-s 


Nv  +  Nv  =  0 


Nxy  +  Njjy  -  0 
Mx  +  Nx®  =  0 


Mx,x  -  (Mxy,y  '  Myx,y)  +  NxW^^  +  NxyW^y  =  Oj 


or 


u  =  0 
V  =  0 
W  Y  =  0 


w  =  0 


(B15) 


The  nonlinear  equilibrium  equations  and  associated  boundary  conditions  (eqs.  (B14) 
and  (B15))  are  used  to  obtain  equations  governing  the  prebuckling  and  buckling  states  of 
the  stiffened  shell  by  taking  the  displacements  u,  v,  and  w  to  be 


u(x,y)  =  ua(x)  +  UB(x,y) 

v(x,y)  =  v^(x)  +  VB(x,y)  >  (B16) 

w(x,y)  =  w^(x)  +  WB(x,y) 

J 

In  equations  (B16)  the  subscript  A  denotes  the  axisymmetric,  prebuckling  displacements 
of  the  stiffened  shell,  and  the  subscript  B  denotes  the  infinitesimal  nonaxisymmetric 
displacements  that  occur  at  buckling. 


Prebuckling  Equations 

For  axisymmetric  prebuckling  deformations,  the  equilibrium  equations  are  found 
from  equations  (B14)  by  discarding  terms  which  have  derivatives  of  y.  The  resulting 
equations  are 
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APPENDIX  B  -  Continued 


^XA,x  =  0 
^-yA,x  =  « 


(B17a) 

(B17b) 


^XA,xx  "W  "  ^XaWa,xx  +  P  -  0 


(B17c) 


with  boundary  conditions 


Nx^  +  Nx  =  0  or 

o 

II 

(B18a) 

Nxy^  +  Nxy  =  0  or 

o 

1! 

<5 

> 

(B18b) 

Mx  +  Nx?  =  0  or 

A. 

wa,x  =  0 

(B18c) 

*A,x  =  » 

or  w^  =  0 

(B18d) 

where  the  subscript  A  terms  are  defined  by  inserting  u  =  Uj^(x),  v  =  v^(x),  and 
w  =  w^(x)  in  the  definitions  of  the  stress  resultants  (eqs,  (B2)  and  (B9)). 

As  discussed  in  the  text,  the  cylinder  was  considered  to  be  free  to  expand  longitudi 
nally  ^u^  #  so  that  integration  of  equation  (B17a)  yields 


Nx^  =  -Nx 


(B19) 


where  Nx  is  the  applied  compressive  load  at  the  boundary.  From  equation  (B19)  and 
the  definition  of  Ny  , 


Nrp  +  Nrp  K-i  ‘t 

- ^+^w 


^11  ^12^A  Nx 

Cii  A,xx-c,,  R  "cii 


(B20) 


The  use  of  the  definitions  of  Mx ,  and  Ny  together  with  equation  (B20)  in  equa- 

A 

tion  (B17c)  yields  the  following  equation  in  w^: 


"i' 


11...  ffff  .  2  ^11^12 


12  A 


-  Kio  +Nv  w”  +^{C 


^  (B21) 
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APPENDIX  B  -  Continued 


Nrr(x)  -  +  Ct, 

Nt(x)-c^^^Ctx  ^2  ^^Ts  ^2  J 


-  kx^eg-/%KT3^e^/-^ 


where  the  primes  denote  total  differentiation  with  respect  to  x. 

Buckling  Equations 

The  equilibrium  equations  and  boundary  equations  which  govern  the  buckling  behaV' 
ior  of  the  stiffened  cylinder  are  obtained  by  substituting  equations  (B16)  into  equa¬ 
tions  (B14)  and  (B15).  K  only  linear  terms  in  the  buckling  displacements  (subscript  B 
terms)  are  retained,  and  if  equations  (B17)  are  utilized,  the  following  buckling  equations 
are  obtained: 


Nxv„  +  Ny  =  0 

^yB,x  yB,y 


yB  - 

Mx  +  Mxv  -  -  M„  +  +  NxWt>  -  Nx»w . 

^B,xx  ^yB,xy  y^B,xy  yB,yy  H  B  A 


(B22) 


-%A'^B,yy  +  ^NxyWg^y  =  0 


with  the  boimdary  conditions 


Nxb  =  0 


Ub  =  0 


NxyB  =  0 


=  0 
^B 


Vb  =  0 


Wb,x  =  0 


(B23) 


^XB,x  ■  (^^B,y  ■  ’^yxB,y 


-  NxWb,x  +  NxbW^  -  NxyWB,y  =  0 


Wb  =  0 


where  the  stress  resultants  with  the  B  subscript  are  given  as  in  equation  (B9)  with  the 
strains  and  curvatures  defined  as 
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^Xb  =  «B,x  +  Vb,x 


yxy^  =  ^  ^B,x  +  VB,y 


'^x  =  ■Wb,xx 


V  -  -'^B,yy 


V  =  -2wB,xy 


(B24) 


It  is  assumed  that  there  are  no  additional  thermal  moments  and  thermal  forces  during 
the  buckling  process  so  that 


Ntxr  =  NTyp.  =  Ntso  =  =  Mtx^  =  MTy  =  M-js^  =  Mxr^  =  0 
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NUMERICAL  SOLUTIONS 

The  solutions  to  the  prebuckling  and  buckling  equations  were  found  by  employing 
the  finite-difference  method.  The  equations  were  formed  into  a  system  of  second-order 
difference  equations.  As  shown  in  reference  25,  this  type  of  system  can  be  easily  solved 
by  matrix  algebra  with  a  modified  Gaussian  elimination  technique. 

In  many  thermal  buckling  problems  of  technical  interest,  a  knowledge  of  the  inter¬ 
action  between  thermal  and  mechanical  loads  is  desired.  Hence,  for  the  present  problem 
there  are  two  possible  choices  for  a  buckling  parameter:  (1)  the  magnitude  of  the  applied 
mechanical  load  (for  example,  the  compressive  load  or  external  pressure)  for  a  fixed 
temperature  distribution  or  (2)  the  magnitude  of  a  temperature  change  which  determines 
thermal  forces  and  moments  for  fixed  mechanical  loading.  In  the  present  study  a  solu¬ 
tion  was  developed  for  each  parameter.  The  solution  using  mechanical  loading  as  the 
buckling  parameter  is  similar  in  form  to  that  of  reference  17  so  that  the  details  are 
omitted  herein.  The  solution  with  temperature  as  the  buckUng  parameter  is  described  in 
detail  in  the  following  sections. 

To  employ  temperature  as  a  buckling  parameter,  a  scalar  temperature  change  T^ 
is  defined  so  that  Tg^  =  '^s  =  and  Tr  =  T^^r;  thus,  in  the  shell 

T  =  TA^sk(l  +  (Cla) 

whereas  in  the  stringers 

T  =  Ta/3s  (l  + 

and  in  the  rings 

T  =  Ta^tI^  + 

Note  that  in  the  integrals  associated  with  equations  (Bll)  and  (B12),  T^  can  be  con¬ 
sidered  as  a  constant. 

By  virtue  of  equations  (Cl),  the  magnitude  of  the  relative  temperatures  in  the  stiff¬ 
ened  cylinder  is  linearly  related  to  a  scalar  T^*  The  shape  of  the  temperature  distri¬ 
bution  is  specified  by  the  constants  /3,  y,  and  g.  Thus,  the  problem  to  be  solved  is  one 
in  which  the  magnitude  of  the  base  temperature  T^  may  vary  while  the  shape  of  the 
temperature  distribution  remains  constant. 
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Prebuckling  Solution 

To  formulate  the  prebuckling  equations  as  a  second-order  finite-difference  system, 
the  variable  ^  is  introduced  so  that 

W;  =  ^  (C2) 


If  the  first  and  second  derivatives  at  the  ith  station  of  the  stiffened  shell  are  approximated 
by  central  differences  as 


(f')i  = 
(f’')i  = 


^i+1  ~  ^i-1 
2A 

ff+l  -  2fi  +  fi_i 


the  prebuckling  equilibrium  equation  (B21)  and  the  definition  of  the  variable  ^  (eq.  (C2)) 
can  be  written  at  the  ith  station  as 

^1  +  ^2^i  +  %(i)(wA)i  =  ^4  +  TAl35(i) 


li  =  0 


^1  =^11  - 


2Fi1^12 


Ry  C 

'11 

^J_ 

r2 

C22( 

Ci2Nx 

-  Ki2  +  Nx 


Cii  R 


1  NTy(i)  +  NTr(i)  -  ^(ntx  +  Nts)  +  ^ 


.gx/a|  -  kr^  ^  +  K.,.  e^^/^ 
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where 

with  the  integration  to  be  performed  over  the  shell  wall  and 

6iH  =1  (i  =  H) 

6iH  =  0  (i  H) 

The  subscript  H  denotes  the  difference  station  at  which  a  ring  is  located.  The  barred 
thermal  terms  are  found  by  considering  a  unit  temperature  change  and  letting 

"^sk  ~  ^sk 

Ts  =  iSg 


Tf  -  ^r 


in  equations  (Bll)  and  (B12). 

In  considering  the  boundary  conditions  associated  with  equation  (C4),  it  was  assumed 
that  the  cylinder  was  supported  at  each  end  (w  =  0)  and  that  the  ends  were  free  to  expand 
longitudinally  and  a  compressive  load  is  prescribed  at  the  boundary.  Under  these 
assumptions 


(C5) 


where  Nx  is  the  magnitude  of  the  applied  compressive  load.  If  the  cylinder  is  supported 
so  that  w^  =  0  at  the  ends,  the  boundary  conditions  on  w  (eqs.  (B18c)  and  (B18d))  can 
be  stated  in  finite-difference  form  at  stations  i  =  0  and  i  =  k  (at  the  ends  of  the  cylin¬ 
der  where  x  =  0  and  x  =  a,  respectively)  as  follows:  For  simply  supported  shells 


^llo  =  ^6  + 

/3l^k  =  %  +  TA^7(kJ 


(C6a) 
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and  for  clamped  shells  (by  using  the  definition  of  |  and  eq.  (C3)) 


where 


and 


-  2(»a)i  '  0 

-  2{wA)k.i  ■ 


^6 


Kll\ 

Cll^ 


^7(l)=^(NTs  +  NTxj- 


M 


Tx 


(C6b) 


(i  =  0,k) 


The  governing  equations  (C4)  and  boundary  conditions  (C6)  can  be  written  in  matrix 
form  as 


L  K 

'  O 

1 _ 

I  Li  I 

(Za)i 

TaRi  +  s 

I  L2  I 

(2^)2 

i 

TaR2  +  S 

I  Lk_i  I 

* 

(ZA)k-l 

^A^k-1  +  ® 

K  L 

1 

_ 1 

TARk  +  S 

where 


(C7) 


(i  =  0,1,2,.  .  .k) 


(i  =  1,2,3,.  .  .k-1) 
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For  simple  support  boundary  conditions 

L  =  I 


and  for  clamped  boundary  conditions 
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Because  equation  (C7)  contains  a  band  matrix,  to  conserve  computer  storage  it  is 
convenient  to  solve  for  the  unknowns  by  Gaussian  elimination.  By  performing  a 

series  of  elementary  row  operations  (see  ref.  17),  equation  (C7)  can  be  written  as 


where 


Pq  =  L-^K 

Pi  =  (Li  -  Pi-l)'* 


(i  =  1.2,. 


(C9) 
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The  unknowns  are  determined  by  finding  Pj^,  Qj,  and  (^Qx)i  from  equa¬ 
tions  (C9),  (CIO),  and  (Cll)  and  solving  successively  for  (z^k-2»  •  •  • 

equations  (C8). 

It  should  be  noted  from  equations  (C8)  that  the  prebuckling  radial  displacements 
are  linear  with  so  that 


+  Ta  ZaJ. 


^i  (%)i 


(C12) 


where  are  the  variables  associated  with  deformation  resulting  from  mechanical 

loading  whereas  ^Z^^^ .  are  the  variables  associated  with  deformation  resulting  from 

thermal  loading.  Deformations  associated  with  mechanical  loading  are  found  by  setting 
equal  to  zero  in  equation  (C8)  and  solving  for  Deformations  associated  with 

thermal  loadings  are  obtained  by  deleting  the  Qj  column  from  the  right-hand  side  of 
equation  (C8)  and  solving  for  corresponding  to  a  unit  base  temperature  change 

(ta  =  l). 

Buckling  Solution 


If  torsional  loading  is  neglected  (N^y  =  Oj,  the  partial  differential  equations  gov¬ 
erning  buckling  of  a  stiffened  cylinder  (eqs.  (B22)  and  (B23))  can  be  uncoupled  into  a  set 
of  ordinary  differential  equations  by  assuming  the  following  relations: 

Ug  =  U(x)  cos  ^  "1 


Vg  =  V(x)  sin 


Wg  =  W(x)  cos 


(C13) 


M(x)  cos 


where  n,  the  number  of  circumferential  buckling  waves,  is  an  integer.  The  substitution 
of  equations  (C13)  into  equations  (B22)  together  with  the  definition  of  Mjjg  yields  four 

ordinary  differential  equations  in  u,  v,  w,  and  M.  To  eliminate  derivatives  of  higher 
order  than  2,  equation  (B20)  was  employed  together  with  the  expression 
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X  ^11  ^12  1  ^12  ^11/  \  ^12 

^B,xxx  ‘“d^  ^  dH  dH  ^  ^  DlirA,x^B,x),x  "  ^Bpcyy 

(C14) 

obtained  from  the  derivative  of  Mx„.  The  resulting  buckling  equations  can  be  written  as 

13 


AiiU"  +  A12U  +  A13V’  +  (A  14  -  TaAjsW"  +  A 16  -  TaA^W' 


+  (Ais  -  T^AigW  +  AiiqM'  =  0 


(Cl  5a) 


A21U’  +  A22V"  +  A23V  +  A24W"  +  (A25  -  TaA26)w'  +  (A27  -  TaA28)w  =  0  (C15b) 


(^31  -  TaA32)u'  +  A33V"  +  (A34  -  TAA35jv  +  (A36  -  TAA37jw’’  +  (A38  -  TAA39 
-  4^310)^’  +  (A311  -  TaA312)w  -  M"  =  0  (C 


(C15c) 


where 


A41U'  +  A42V  +  A43W"  +  ^A44  -  TaA45)w'  +  A46W  +  M  =  0  (C15d) 


All  =  Cii 


A 12  =  -^66  n 


Ai3  =  ^12  +  C66 


Kl2Kiin 


Dll  R 


Ai4  =  AllwV 


Ai5  =  -Aiiw]^^ 


A  _  A  ,  1 L  ^11^12)  ,  L-  01^  KiiDi2|/n\2 

+  bR2  -  +  Ki2  ♦  2K66  -  -5— ) fej 
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A17  = 

^18  =  -‘^66(|j  ^Am 

^19  =  ^At 

"^110 '^n 

A2i  =  -(Ci2  +  Cee)^ 
A22  =  ^66 

^  /n\2 

A23  =  “^22fe 

A24  =  ^Kj2  +  2K0^^ 
A25  =  ^21^A]^ 

A26  =  -A21^iT 


A27  =  -C66  % 


C22^ 
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A41  =  -Kii 

A42  =  -Ki2  I 

^43  =°11 
A44  = 

A45  =  Kiiw^^ 


''^46 


To  represent  the  Dirac  delta  and  doublet  functions  appearing  in  the  coefficients  of 
equation  (CIS),  the  following  finite-difference  approximations  were  employed: 

6(x  -  ji)  f  =  6i|j  ^ 


6'(x  -jl)t  =  6(i-i)H^  -  6(i+i)H^ 

In  this  form,  equations  (CIS)  are  the  analogs  of  equations  (36)  of  reference  17,  As 
outlined  in  that  reference,  equations  of  this  form  can  be  cast  into  matrix  form,  and  buck¬ 
ling  loads  can  be  found  by  examination  of  a  simple  determinant  resulting  from  a  Gaussian 
elimination  procedure  similar  to  that  previously  described  herein.  A  solution  of  this  type 
was  obtained  in  the  present  study  by  using  the  determinant  plotting  method  described  in 
reference  17. 


Modal  Iteration  Scheme 

If  the  central  difference  formula  (eqs.  (C3))  are  employed  for  the  first  and  second 
derivatives  of  U,  V,  W,  and  M  in  equations  (CIS),  the  equilibrium  equations  can  be 
written  at  the  interior  stations  of  the  shell  as 


AiZi_i  +  BiZj  +  CjZj^j  =  T^^D^Zj,!  +  E^Zj  +  FiZ^^j 

+  T2(GiZi.i  +  HiZi^.i 


(i  =  1,2,.  .  .k-1)  (C16) 
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where 


All 

-fAi3 

Ai4  "f  Ai6 

-f  Alio 

-f  A21 

A22 

A24  -  J  A25 

0 

-|a31 

A33 

A36  -f  A38 

-1 

-f  A41 

0 

A43  -  f  A44 

0 

-2Aii  +  a2Ai2 

0 

-2Ai4  +  A^Ais 

0 

B4  = 

0 

-2A22  +  a2a23 

-2A24  +  a2A27 

0 

1 

0 

-2A33  +  a2a34 

-2A30  +  a2a3ii 

2 

0 

A^A42 

-2A43  +  A2A4e 

a2 

All 

f  Ai3 

Ai4  A10 

A  A 

^Ano 

f  A21 

A22 

A24  +f  A25 

0 

fA31 

A33 

A36  +  f  A38 

- 

■1 

|a41 

0 

A43  +  f  A44 

0 

0 

0 

A15  ■f’ Ai7 

0 

0 

0 

-f  A26 

0 

D.  = 

1 

-f  A32 

0 

A37  -  f  A39 

0 

0 

0 

-f  A45 

0 
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In  order  to  develop  a  modal  iteration  solution,  the  following  forward  and  backward 
difference  formulas  are  introduced  at  ends  of  the  cylinder  (finite-difference  stations  0 
and  k) 


(C17) 


The  cylinder  was  assumed  to  be  supported  at  the  ends  so  that  Wg  =  0.  The  eight  possible 
boundary  conditions  for  buckling  in  equations  (B23)  can  be  written  in  matrix  form  as 


BqZo  +  CqZi  =  T^EqZj 
•^k^'k-l  +  ®k^k  =  "^A^k^k-l 


(CIS) 


with 


-All 

AnL  KiiKi2\ 

j 

0 

Dll 

1  -  0-11 

0 

0 

0 

-C66|a 

-^66 

0 

0 

+ 

0 

1  -'^2 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-1 

0 

0 

0 

0 

1 

pi 

_ 

(C19a) 


0 


^ll^A 


M 


(0) 


0 


0 


^66 

0 


2K66- 


R 


0 


0 

0 


(C19b) 


0 


0 


1 


0 
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where 


and 


Eq  =  a 


0  0 


0 


a  = 


0 

0 


0  0 


0 

0 


0 


0  0 


“44 


Ak  =  Co  1 


Bk  =  Bq  ) 


Ek  =  -EqJ 


(C19c] 


(C20) 


where  in  equations  (C20)  w^  terms  are  to  be  evaluated  at  x  =  a  rather  than  at  x  =  0. 
The  elements  of  the  selection  matrix  a  are  defined  as  follows: 


Element 

Value 

Boundary  condition 

“11 

0 

Ub  =  0 

1 

Nxt,  =  0 

Xb 

^2 

0 

Vb  =  0 

1 

“44 

0 

Mx  =0 

^B 

1 

Vb  =  0 

APPENDIX  C  -  Concluded 


The  equilibrium  equations  (C16)  and  boundary  conditions  (C18)  can  now  be  written 
in  matrix  form  as 

AX  =  T^BX  +  t|cX  (C21) 

where  is  a  scalar,  X  is  a  vector,  and  A,  B,  and  C  are  matrices.  In  order  to 
perform  modal  iteration,  let 

(C22) 

Ta 

so  that  equation  (C21)  can  be  written  as 

XAX  =  BX  +  1  CX  (C23) 

To  start  the  iteration  process  assume  that  X  =  X°  where  X*^  is  a  trial  eigenvector. 
Substitute  X®  into  the  right-hand  side  of  the  approximate  equation 

AAX  *  BX 

and  solve 

AAX  =  BX°  (C24) 

by  treating  AX  as  the  unknown  variable.  In  finding  the  solution  to  equation  (C24) , 
Gaussian  elimination  (Potter's  method)  is  employed  to  avoid  having  to  find  the  inverse 
of  A.  The  resulting  approximate  eigenvalue  Aj  and  normalized  eigenvector  X^  can 
be  taken  as 

Ai=lAXl 

Xl  =  M  (C25) 

For  subsequent  iterations  the  equation 

AAX  =  BXi  +  —  CXi  (i  =  1,2,.  .  .)  (C26) 

must  be  solved  to  find  more  accurate  solutions  to  equation  (C23),  The  iteration  process 
converges  to  the  largest  characteristic  value  of  X  (lowest  T^)  provided  —  CX  is 

small  compared  with  BX.  Thus  the  temperature  T^  can  be  found  for  a  specified  num¬ 
ber  of  circumferential  waves  n.  The  lowest  value  of  T^  found  is  then  the  buckling 
temperature. 
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COMPUTER  PROGRAMS 
BAMSOC  I 

The  computer  program  BAMSOC  I  finds  the  buckling  load  for  a  stiffened  cylindrical 
shell  subjected  to  mechanical  and  thermal  loadings  by  using  determinant  plotting.  The 
programed  solution  is  similar  to  that  developed  in  reference  17  and  differs  only  in  that 
the  present  solution  considers  thermal  effects  in  a  multilayered  orthotropic  cylinder. 
Buckling  loads  or  buckling  temperatures  can  be  selected  as  the  buckling  parameter,  and 
critical  values  are  found  by  determinant  plotting  over  user -specified  ranges  of  circum¬ 
ferential  wave  number  n  and  the  buckling  parameter.  For  a  specified  n,  the  program 
seeks  a  change  in  sign  in  a  determinant  or  "modified  residual"  by  varying  the  buckling 
parameter.  A  flow  diagram  for  BAMSOC  I  (determinant  plotting  solution)  follows: 
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The  program  is  written  in  FORTRAN  IV  and  will  run  on  the  Control  Data  6000 
series  computer  with  the  Scope  3.0  operating  system.  The  program  elements  were  over- 
layed  so  that  a  maximum  of  200  finite-difference  intervals  could  be  used  on  a  remote 
computer  terminal  with  70  000  octal  storage  units.  A  maximum  of  20  layers  can  be  used 
in  the  shell.  Input  required  by  the  program  is  contained  in  the  following  namelists: 

$  CON 

NI  initial  value  of  n  to  be  investigated 


NF 

DELN 

NBARXI 

NBARXF 

DELNX 

PINIT 

PFIN 

DEEP 

TAI 
TAF 
DELTA 
EASE  =  1 
=  2 

=  3 
=  4 


final  value  of  n 
increment  in  n 

initial  or  fixed  value  of  N^,  axial  compressive  load 
final  value  of 
increment  in 

initial  or  fixed  value  of  p  or  q,  lateral  or  hydrostatic  pressure 
final  value  of  p  or  q 
increment  in  p  or  q 

initial  or  fixed  value  of  base  temperature  (see  eq.  (Cl)) 
final  value  of 
increment  in 

axial  compression  with  fixed  temperature  and  lateral  pressure 
lateral  pressure  with  fixed  temperature  and  axial  compression 
hydrostatic  pressure  with  fixed  temperature  and  axial  compression 
thermal  loads  corresponding  to  with  fixed  mechanical  loads 


$  (end  of  CON) 
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$  NANCY 
R 

A 

DELBAR 
BC  =  1 

=  2 
ER 
AR 

IR 

ZR 

GRJR 

NRING 

NRSPACE 

ES 

AS 

IS 

zs 

GSJS 

D 


radius  to  reference  surface 
cylinder  length 

distance  from  inner  wall  surface  to  reference  surface 
simply  supported  for  prebuckling  analysis 
clamped  for  prebuckling  analysis 
Young's  modulus  for  ring 
area  of  ring 

moment  of  inertia  of  ring  about  its  centroid 

distance  of  centroid  of  ring  to  reference  surface  (positive,  external  rings) 
ring  torsional  stiffness,  GJ 
number  of  rings 

number  of  finite-difference  intervals  between  adjacent  rings  (no  more 
than  200  stations  along  cylinder  length) 

Young's  modulus  for  stringer 

area  of  stringer 

moment  of  inertia  of  stringer  about  its  centroid 

distance  from  centroid  of  stringer  to  reference  surface  (positive, 
external  stringers) 

stringer  torsional  stiffness,  GJ 

stringer  spacing 
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ECC  distance  from  reference  surface  to  surface  on  which  Nx  ^^ts 

GAMS  stringer  temperature  shape  parameter,  (see  eq.  (Clb)) 

GAME  ring  temperature  shape  parameter,  y^  (see  eq.  (Clc)) 

TAUS  J  0!g^s(z)dAs 

TAUR  y  Q!j.^j,(z)dAr  where  /3  is  temperature  shape  parameter  in 

„  )  equation  (Cl) 

TAUSPR  \  Q!s|3s(z)z  dAg 

TAURPR  r  a^Pj.{z)z  dA^ 

d  J 

G  exponential  temperature  shape  factor,  g  (see  eq.  (Cl)) 

EX(I)  Ex  of  ith  layer  of  shell  (first  layer  lies  on  inner  surface  of  shell;  no 

more  than  20  layers) 

EY(I)  Ey 

NUX(I)  Px 

NUY(I)  My 

GXYtt)  Gxy 

H(I)  thickness  of  layer 

GAMMA(I)  layer  temperature  shape  factor  ygj^  (see  eq.  (Cla)) 

BETA(I)  layer  temperature  shape  factor  /3sk  (Cla)) 

ALPHAX(I)  ax 

ALPHAY(I)  Oy 

LAYER  number  of  layers  (=20) 

$  (end  of  NANCY) 
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$  CORNIE 

lALF  =  1 

Nxb  =  0, 

Vb  =  0,  Wb  =  0,  Mxg  =  0 

=  2 

Ub  =  0, 

Vb  =  0,  Wb  =  0,  Mxg  =  0 

=  3 

Nx„  =  0, 

Nxyg  =  0,  Wb  =  0,  Mxg  =  0 

=  4 

Ub  =  0, 

NxyB  =  0,  Wb  =  0,  Mxg  =  0 

=  5 

Ub  =  0, 

^xy^-0,  Wb  =  0,  Wb^x  =  0 

=  6 

Nxb  =  0, 

NxyB  =  0’  '^B  =  0»  ^B,x  =  ® 

=  7 

o' 

11 

PD 

Vb  =  0,  Wb  =  0,  Wbx  =  0 

=  8 

Ub  =  0, 

Vb  =  0,  Wb  =  0,  Wbx  =  0 

MODE  =  0  calculation  of  buckling  mode  shape  is  not  required 

4-  0  calculation  of  buckling  mode  shape  is  required 

IRING  =  1  simplified  ring  theory  described  in  appendixes  A,  B,  and  C 
=  2  more  accurate  ring  theory  described  in  appendix  E 
$  (end  of  CORNIE) 

$  RING  (required  only  if  IRING  =  2) 

IX  Ix,  moment  of  inertia  for  out-of-plane  bending  of  ring  (eq.  (Ell)) 

lOXZ  loxz)  cross  product  of  inertia  for  out-of-plane  bending  of  ring  (eq.  (Ell)) 

$  (end  of  RING) 

A  sample  output  listing  for  BAMSOC  I  follows.  The  case  shown  is  the  ring-  and 
stringer-stiffened  cylinder  loaded  by  thermal  loading  alone.  (See  text.)  The  parameter 
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has  been  selected  as  the  change  in  temperature  of  the  heated  skin.  Experience  with 
determinant  plotting  solutions  has  shown  that  it  is  possible  to  miss  sign  changes  in  the 
determinant  by  taking  too  large  an  increment  in  the  load  parameter.  As  a  safeguard 
against  this  possibility,  the  present  program  examines  small  loading  increments,  when 
unusual  trends  in  the  determinant  are  detected. 

In  the  sample  problem,  for  N  =  51  or  52,  the  modified  residual  value  was  found  to 
increase  in  value  after  a  steady  decline.  The  irregularity  was  detected  and  finer  load 
increments  were  employed  in  the  calculations  and  a  change  in  the  sign  of  the  determinant 
was  found.  Note,  however,  that  no  irregularities  were  found  for  N  =  53,  and  as  a  con¬ 
sequence,  sign  changes  were  found  at  higher  temperatures.  Further  investigation  with 
finer  load  increments  in  calculations  not  presented  established  a  buckling  temperature  of 
about  525°  F  (292  K)  for  N  =  57. 
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namelist  input  for  BAMSOC  1 


SCON  NIsSO,,  NF«53.«  DELN=1.»  NBARX!=0*«  NBARXF=0*»  DELNX=1.. 
PIN1T=0.*  PFIN=0*«  DELP=1.«  TA1=450.«  TAF=600»i  DELTA=20«*  KASE=4S 


SNANCY  R=65.,  A=300.  ,  DELBaR  =  .03,  BC=1.  ER  =  16.4E6.  ‘  ^ 

ZR*-I.992.  GRJR=3477.»  NRING=19,  NRSPACE=10.  ES=16.4E6»  AS=.C495.  IS-. 004082. 

2S=“.342*  GSJS=93..  D=1.35.  ECC=0..  GAMS=0. «  GAMR-0..  TAUS=0..  TAUR=0.' 
TAUSPR  =  0..  TAURPR  =  0..  G  =  0.,  EX  (  1  )  =  ’  4 . 5E6 .  EY  (  1  )  =  1  4 . 5E6 ,  NUX  (  1  )  =  .  3  .  NUY(l)-.->. 
GXY( 1 )=5»580E6«  H(l)=.06.  GAMMA(1)=0..  BET A ( 1 ) = 1 . . 


ALPHAXd  )=*5.E-6.  ALPHAY  ( 1  )  =5 . E-6  *  LAYER=1$ 
SCORNIE  IALF=1*  M0DE*0.  iRfNG^lE 


BAMSOC  I 

CALCULATION  OF  BUCKLING  LOAD  FOR  THERMALLY  STRESSED,  STIFFENED,  MULTI-LAYERED  CYLINDER  USING  DETERMINANT  PLOTTING 
CHANG-CARO  R0F364  NASA  LANGLEY  RESEARCH  CENTER  1968 

DATE  03/27/69 

NI*  5.0000E+01  NF=  5.3000E+01  DELN=  l,0000E+00 

TAI-  4.5000E+02  TAF=  6.0000E+02  DELTA=  2.0000E+01 

NBARX=  0. 

P=  0. 


INPUT  FOR  PROGRAM  THERMUL 


CHANG-CARD  A2031  RDF364  JANUARY, 1968 


GEOMETRY  PARAMETERS  LOADINGS 

R=  6.5000E^-01 
A=  3.0000E<-02 

DELBAR*  3.0000E-02  ECC=  0. 


TEMPERATURE  PARAMETERS 

TAUS=  0.  GAMS=  0. 

TAUR=  0.  GAMR=  0. 

TAUSPR=  0.  G=  0. 

TAURPR=  0. 


RING  PROPERTIES 

ER=  1.6400E+07  ZR=-1 .9920E+00 
AR=  2.4600E-01  GRJR=  3.4770E+03 
IR=  1.3770E-01 


STRINGER  PROPERTIES 
ES=  1.6400E<-07 
AS=  4.9500E-02 
IS=  4.0820E-03 


ZS=-3.4200E-01 
GSJS=  9.3000E+01 
0=  1.3500E+00 


NUMBER  OF  RINGS=  19 

NUMBER  OF  FINITE  DIFFERENCE  SPACES  BETWEEN  RINGS=  10 


GEOMETRY  AND  TEMPERATURE  PROPERTIES  FOR  LAYERS 


LAYER 

1 


EX  EY  NUX 

1.450E+07  1.450E+07  3.000E-01 


NUY 

3.000E-01 


GXY  H  GAMMA 

5.580E+06  6,OOOE-02  0. 


BETA  ALPHA  X  ALPHA  Y 
l.OOOE+00  5.000E-06  5.000E-06 


SIMPLE  SUPPORT  BOUNDARY  CONDITIONS 


11 


12 


22 


66 


13 


23 


1.59340659E+07 
9.56043956E+05 
2.86813187E+02 

0. 

L 

M 


5.58000000E+06  -1 .03571429E+02  -1 • 0357 1429E+02 

3.34800000E+05 

1,00440000E+02 

0. 

-6,21428571E+00  -6.21 428571E+00 

0.  0. 


B  1  1.593406596+07  4.7802 1978E+06 
C  9.56043956E+05  2. 86813187E+05 
D  2.86813187E+02  8.60439560E+01 
K  0.  0. 
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INPUT  FOR  CALCULATION  OF  MODIFIED  RESIDUAL  CHANG-CARO  R0F364  FEBRUARY, 1 968 

IALF=  1 


N=  5.C000E-I-01 

TA  MODIFIED  RESIDUAL 

4.50000000E+02  -7 .79225885E+05 

4.70000000E+02  -7 .75252974E *05 

4.90000000E+02  -7 . 65576861E+05 

5.10000000E+02  -7 •60571046E*05 

5.30000000E4-02  -7  •42188514E+05 

5.50000000E+02  6 .7900992lEf 05 

5.40000000E+02  -7 . 14898570Ef 05 

5.40000000E+02  CROSSING  FOUND  WITHIN  1.  PERCENT 


N=  5.1000E+01 

TA  MODIFIED  RESIDUAL 


4.50000000E+02 
4.70000000E+02 
4.90000000E+02 
5.10000000E+02 
5.30000000E+02 
5.50000000E+02 
OPEN  UP  INTERVAL  IN  SEAR 
5.30000000E+02 
5.32000000E+02 
5.34COOOO0E+O2 
5.36000000E+02 
5.38000000E+02 
5.40COOOOOE-I-02 
5.39000000E+02 
5.39000000E+02 


-7.88803726E+05 
-7.844748i9E+05 
-7.78172671E+05 
-7.67786896E«-05 
~7.43541985E+05 
-9.24281493E+05 
CH  OF  CROSSING 

“7.43541985E+05 
-7.38292392E<'05 
“7.31168334Ei-05 
-7.19848153E«-05 
“6.9I640017E+05 
9.03024088E+05 
-6.24795117E+05 
CROSSING  FOUND  WITHIN 


PERCENT 


5.20COE  +  01 


TA 


MODIFIED  RESIDUAL 


4.50000000E+02 
4. 70000000 E+02 
4.90000000E+02 
5.10Q0O000E+02 
5.30000000E+02 
5.50000000E+02 
OPEN  UP  INTERVAL  IN  S 
5.30C000O0E+02 
5.32000000E+02 
5.34000000E+02 
5.36000000E+02 
5.35000000E+02 
5.35000000E+02 


-7.98255606E+05 
~7,93594324E*-05 
^7.86682408£+05 
-7.74838  733E*-0  5 
-7.41093951E+05 
-7.61622322E+05 
CH  OF  CROSSING 

-7.41093951E+05 
-7.2g620508E+05 
-6.99861698E+05 
8.61174311E+05 
“6.16100274E+05 
CROSSING  FOUND  WITHIN 


PERCENT 


N=  5.3000E+01 

TA  MODIFIED  RESIDUAL 

4.50000000E-I-02  -8 . 07594308E«-05 

4.70000000E-I-02  -8 .02630842E  +  05 

4.9C000000E+02  -7  .95147042E4-0  5 

5.10C00000E+02  -7.81817375E+05 

5.30000000E+02  -7 .24361 857E+05 

5 .50000000E+02  -7 ,05664617Ef05 

5.70000000E  +  02  6 .79333836E<-05 

5.60000000E+02  8 .50808  776E«-05 

5.60000000E+02  CROSSING  FOUND  WITHIN  1.  PERCENT 


MINIMUM  TA=  5.35000000E+02  N=  5 . 20000000E+0 1 
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BAMSOC  n 

The  computer  program  BAMSOC  II  finds  the  buckling  load  for  a  stiffened  cylindrical 
shell  subjected  to  mechanical  and  thermal  loadings  by  using  modal  iteration.  The  solu¬ 
tion  corresponds  to  equations  (C16),  (C18),  and  (C27)  of  the  present  text.  A  flow  diagram 
for  BAMSOC  11  (modal  iteration  solution)  follows: 


Buckling  temperatures  are  found  for  specified  ranges  of  the  circumferential  wave  num¬ 
ber  n.  If  the  solution  is  well  behaved,  the  program  will  converge  to  the  candidate  buck¬ 
ling  temperature  after  several  iterations.  The  program  is  overlayed  to  permit  a  maxi¬ 
mum  of  200  finite-difference  intervals  and  requires  150  000  octal  storage  units. 
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The  input  required  by  the  modal  iteration  solution  is  similar  to  that  for  BAMSOC  I 
except  that  loading  ranges  and  increments  are  not  specified  and  the  solution  is  valid  only 
for  simply  supported  boundary  conditions.  The  following  modifications  are  required  for 
the  three  namelists  presented  for  BAMSOC  I: 

(1)  in  namelist  CON,  retain  only  NI,  NF,  DELN,  NBARXI,  and  PINIT 

(2)  retain  namelist  NANCY 

(3)  in  namelist  CORNIE,  permissible  values  of  lALF  are  from  1  to  4.  Delete 

MODE  and  IRING. 

The  convergence  criteria  (see  eq.  (C27))  is  specified  internally  for  the  ith  iteration  as 


IW 

i  -  (t. 

A)i-l| 

(Ta)i 

0.001 


A  sample  output  listing  for  BAMSOC  n  follows.  The  case  is  the  same  case  as  for 
BAMSOC  I  except  that  only  100  finite-difference  intervals  were  employed  to  condense  the 
output. 
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NAMtLiST  INPUT  FOR  BAM50C  II 


SCON  NI=57,,  NFsS7,»  DELN=1.*  NdARXI=0.,  PlNlT=0.S 

SNANCY  P=6b.f  As30().*  Qt.L8AR  =  .<jJ,  dC=lf  ER=:16.4E6f  AR=.246»  IR=.1377» 
ZR=-1.992»  GkJR=3477..  NRINO=iy»  NRSPACE=  5»  ES=16.4E6,  AS=.0495,  IS=,004082, 
ZS=-,342,  GSJS=^3.*  0=1.35*  LCC=0.»  GAmS=0.*  GAMR=0,*  TAUS=0.*  TAuR=0,* 
TAUSPR=f^..  TAUROR^O.*  G=0,*  GXY(l)=b.5b0t6,  H(l)  =  .d6,  GAMK1A  ( 1 )  =0  .  *  aETA<l)=l.* 
EX(1)=14.5E6,  EY<l)=14.bE6*  NuA(l>=.30*  NUY(l)  =  .3iJ,  ALPHAX  ( I )  =5  ,E-6, 

ALPHAY ( 1 ) =b.t-6*  LAYEk=1S 

SCORNIE  IAi.F  =  l<s 


BAMSOC  II 

CALCULATION  OF  BUCKLING  LOAD  FOR  THERMALLY  STRESSED,  STIFFENED,  MULTI-LAYER  CYLINDER  BY  MODAL  ITERATION 


C HANG-CARD 

RDF364 

NASA  LANGLEY 

RESEARCH  CENTER 

1968 

DATE  03/27/69 

NI=  5.7000E+01 

NF= 

5.7C00E+01 

DELN=  l.OOOOE+00 

INITIAL  Z  VECTOR  (UNSCALEDI 


UB  VB 

I  0,  0. 

21  0*  0. 

41  0,  0. 

61  0.  0. 

81  0.  0. 

101  0.  0. 

ITERATION  NUMBER  I 


TEMPERATURE  TA=  1 , 63659058E-C 1 


W8  MB 

l.OOOOOOOOE^-OO  0. 
l.OOOOOOOOEt-00  0. 
l.OOOOOOOOE+00  0. 
l,00000000E<-00  0. 
l.OOOOOOOOE+00  0. 
l.OOOOOOOOE+00  0. 


CORRESPONDING  SCALED  Z  VECTOR 


UB 

VB 

MB 

MB 

1 

-7.76570953E-07 

0. 

0. 

0. 

21 

-1.91i30069E-10 

1.83146093E-08 

7.35752664E-09 

-1.83573040E-01 

41 

-7.01489655E-14 

l,83760317E-08 

7.39760844E-09 

-1.83772776E-01 

61 

7.01489688E-14 

1.83 760317E-08 

7.39760844E-09 

-1.83772776E-01 

81 

1.91130C69E-10 

1.83146093E-08 

7.35752664E-09 

-1.83573040E-01 

101 

7.7657C953E-07 

0. 

0. 

0. 

ITERATION  NUMBER  2 


TEMPERATURE  TA=  1 . 17890697E+0 3  CORRESPONDING  SCALED  Z  VECT3R 


UB 

VB 

MB 

MB 

1 

7.27927938E-07 

0. 

0. 

0. 

21 

-1.57725826E-09 

5.04390724E-08 

3.52955211E-08 

-1.52729564E-01 

41 

-1.20398426E-12 

5.07973257E-08 

3.55277678E-08 

-1. 53943644 E-01 

61 

1-20398423E-12 

5.07973257E-08 

3.55277678E-08 

-1.539436446-01 

81 

1.57725826E-09 

5.04390724E-08 

3.52955211E-08 

-1.52729564E-01 

101 

-7-27927938E-07 

0. 

0. 

0. 

ITERATION  NUMBER  25 


TEMPERATURE  TA=  5.37719505E4-02  CORRESPONDING  SCALED  Z  VECTOR 


UB 

VB 

MB 

MB 

1 

4.58064816E-07 

0. 

0. 

0. 

21 

-2.876251696-07 

-9.81722985E-09 

-6.56022956E-09 

3.23334452E-02 

41 

-4. 84982682 E-08 

-4.40150651E-09 

-2.92732002E-09 

1.45714341E-02 

61 

4.84982682E-08 

-4.40150651E-09 

-2.92732002E-09 

1.45714341E-02 

81 

2.87625169E-07 

-9.81722985E-09 

-6.56022956E-09 

3.23334452E-02 

101 

-4.58064816E-07 

0. 

0. 

0. 
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ITERATION  NUMBER  26 


TEMPERATURE  TA=  5 . 37123003E+02  CORRESPONDINS  SCALED  Z  VECTOR 


UB 

VB 

WB 

MB 

I 

4.46578687E-07 

0. 

0. 

0. 

21 

-2.930478268-07 

-9.38196491E-09 

-6.26962170E-09 

3. 08988773 E-02 

41 

-5.299828578-08 

-4.69946566E-09 

-3.126I3561E-09 

1.55542013E-02 

61 

5.299828576-08 

-4.69946566E-09 

-3.12613561E-09 

1.555420136-02 

81 

2.93047826E-07 

-9.38196491E-09 

-6.26962170E-09 

3.08988773E-02 

101 

-4.46578687E-07 

0. 

0, 

0. 

ITERATION  NUMBER  27 


TEMPERATURE  TA=  5.36575924E4-02  CORRESPONDINS  SCALED  Z  VECTOR 


UB 

VB 

WB 

MB 

1 

4.35668C83E-07 

0. 

0, 

0. 

21 

-2.97920426E-07 

-8,949l0655E-09 

-5.980573946-09 

2.94725094E-02 

41 

-5.75I10986E-08 

-4.98896109E-09 

-3.319328516-09 

1.55089071E-02 

61 

5.75110986E-08 

-4.98896109E-09 

-3.31932851E-09 

1.65O09O71E-O2 

81 

2.97920426E-07 

-8.94910655E-09 

-5. 980573 94E-09 

2,947250946-02 

101 

-4.35668C83E-07 

0. 

0. 

0, 

ITERATION  NUMBER  28 


CONVERGED  TEMPERATURE  TA=  5. 3607A232E+02  CORRESPONDING  SCALED  Z  VECTOR 
UB  VB  WB  MB 


1 

4.25301149E-07 

0. 

0. 

2 

-7.59850632E-07 

-4. 59625262 E-07 

2.88703579E-05 

3 

-5.21414819E-07 

-7.71134756E-07 

4.27008188E-05 

4 

7.46228255E-11 

-7,28756070E-07 

3.92369407E-05 

5 

2.70353608E-07 

-4. 414841 70E-07 

2.312306666-05 

6 

2.99574280E-07 

-9.37078924E-09 

-5. 52645553E-09 

7 

2,81086829t-07 

4,287266276-07 

-2,328198626-05 

8 

1,117223316-07 

6.970146706-07 

-3.796409956-05 

9 

-1.45090817E-07 

7.09547922E-07 

-3.85226372E-05 

10 

-3.27532512E-07 

4,551464456-07 

-2.42795296E-05 

11 

-3.62973005E-07 

2.43286648E-09 

1.72226769E-09 

12 

-3.44911857E-07 

-4.67981720E-07 

2. 505843 71E-05 

13 

-1.44970992E-07 

-7.465562186-07 

4.04388090E-05 

14 

1.49747724E-07 

-7.43563747E-07 

4.02956732E-05 

15 

3.451138536-07 

-4.61421596E-07 

2,478434876-05 

16 

3.53028614E-07 

4.79054708E-09 

3.21738900E-09 

17 

3. 03512382 E-07 

4.42124012E-07 

-2.349386756-05 

18 

1. 11311952 E-07 

6. 841915436-07 

-3.705585466-05 

19 

-1.52077748 E-07 

6.67749535E-07 

-3.624634556-05 

20 

-3,162784126-07 

4.04730815E-07 

-2.18899554E-05 

21 

-3.02289280E-07 

-8.52149244E-09 

-5.69499266E-09 

22 

-2.40603980E-07 

-3.70626319E-07 

1.958333956-05 

23 

-7.79210795E-08 

-5. 60691973 E-07 

3.036334536-05 

24 

1.33393533E-07 

-5.38503291E-07 

2.927025756-05 

25 

2.5071O335E-O7 

-3.20106295E-07 

1.74149936E-05 

26 

2.34182913E-07 

9.58375632E-09 

6.397016036-09 

27 

1.73529469E-07 

2.81843089E-07 

-1.481794666-05 

28 

4.90093406E-08 

4.17721507E-07 

-2.262101936-05 

29 

-1.05031425E-C7 

3.95379757E-07 

-2.151999606-05 

30 

-1.92160659E-07 

2.30933170E-07 

-1.263219876-05 

31 

-1.65745567E-07 

-8.73802977E-09 

-5.826604316-09 

32 

-1.14405427E-07 

-1.96104461E-07 

1.02615469E-05 

33 

-2.73157633E-08 

-2.8492  8388E-07 

1.54304282E-05 

34 

7.54101910E-08 

-2.65806306E-07 

1.448774806-05 

35 

1.30679669E-07 

-1.52498812E-07 

8.38802829E-O6 

36 

1.07290082E-07 

7.01425061E-09 

4,672488456-09 

37 

6.81917701E-08 

1.24588002E-07 

-6.48369788E-06 

38 

1.24C77770E-08 

1.76889287E-07 

-9.57024365E-O6 

39 

-4.98547675E-08 

1.62C46387E-07 

-8.84625794E-06 

40 

-8.10843369E-08 

9.07I49850E-08 

-5.02850473E-06 

0. 

1.97151230E-01 
2.03472475E-01 
1.45920054E-*01 
6.134e3855E~02 
‘6.63765673E-04 
-1,08527967E-01 
-1.76416032E-01 
-1.84431761E-0I 
-1.24870198E-01 
-9.63841827E-03 
1.20262254E-0I 
1.9Z4736I3E-01 
1.90551988E-01 
1.15363571E-01 
-1.57017051E-02 
-1.23089686E-01 
-1,78314502E~O1 
-1.68109183E-01 
-9.40429520E-02 
2.80636251E-02 
1.08943509E-01 
1.47419054E-01 
1.33661345E-01 
6.96391649E-02 
-3.15968803E-02 
-8.67100390E-02 
-1.10635773E-01 
*»*9. 67952493  E-02 
-4.70731332E-02 
2.88428069E-02 
6,28887408E-02 
7.59949709E-02 
6.41582849E-02 
2.89144603E-02 
-2.31792101E-02 
-4.18076618E-02 
“4,761483966-02 
-3.04342699E-O2 
-1.53969322E-02 
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41 

-6.20122302E-08 

-5.26917618E-09 

-3,50634900E-09 

1.74328862E-U2 

42 

-3.40867935E-08 

-6.97854476E-08 

3.59568563E-06 

2.52369502E-02 

43 

-2.04212480E-09 

-9,50531674E-08 

5.14159204E-06 

2.61273085E-O2 

44 

3.04004734E>08 

-8.3  8477792  E-*08 

4.58879044E-06 

1.92021907E-02 

45 

4.40193554E-08 

-4.41941266E-08 

2.49591889E-06 

5.26846401E-03 

46 

2,79056634E-08 

4.04701034E-09 

2.69043138E-09 

-1.34039223E-02 

47 

7,88264544E-09 

2.8146t217E-08 

-1.39425150E-06 

-1.29366908E-02 

48 

-6,66980972E-09 

3.22525792E-08 

-1.73154237E-06 

-9.83661853E-03 

49 

-1.66255a69E-08 

2.30691997E-08 

-1.27836897E-06 

-4.16443057E-03 

50 

-1.61454443E-08 

7.16G82486E-09 

-4.92203515E-07 

3.44703134E-03 

51 

-5.11607606E-20 

-3.61446716E-09 

-2.40179628E-09 

1.19772649E-02 

52 

1.61454443E-08 

7.16082466E-09 

-4.92203515E-07 

3.44703134E-03 

53 

1.66255e69E-08 

2.30691997E-08 

-1.27836897E-06 

-4.16443057E-03 

54 

6.66S80972E-09 

3.22525792E-08 

-1.73154237E-06 

-9.83661853E-03 

55 

“7.88264544E-09 

2.81461217E-08 

-1.39425150E-06 

-1.293669O0E-O2 

56 

-2.79055634E-08 

4.04701034E-09 

2.&9043138E-09 

-1.34039223E-02 

57 

-4.40193554E-08 

-4.41941266E-08 

2.49591889E-06 

5.26846401E-03 

58 

-3.04004734E-08 

-8.38477792E-08 

4. 58879044E-06 

1.92021907E-02 

59 

2.04212460E-09 

-9.'50531674E-08 

5.14159204E-06 

2.61273885E-02 

60 

3.40867935E-08 

-6.97854476E-08 

3.595665638-06 

2.52369502E-02 

61 

6.20122302E-08 

-5.26917618E-09 

-3.50634900E-09 

1. 74328862 E-02 

62 

8.10843369E-08 

9.07149850E-08 

-5.02850473E-06 

-1.53969322E-02 

63 

4.98547675E-08 

1.62046387E-07 

-8.84625794E-06 

-3.84342699E-02 

64 

-1. 2407777 OE-08 

1.76889287E-07 

-9.57824365E-06 

-4.76148396E-02 

65 

-6.81917701E-08 

1.24580OO2E-O7 

-6.48369788E-06 

-4.18076618E-02 

66 

-1,07290C82E-»07 

7.01425061E-09 

4.67240845E-O9 

-2.31792101E-02 

67 

-1.30679669E-07 

-1.52498812E-07 

8.388828298-06 

2.89144603E-02 

68 

-7.54101910E-08 

-2.658063C6E-07 

1. 44877480 E-05 

6.41582849E-02 

69 

2.73157633E-08 

-2.84928388E-07 

1.54304282E-05 

7.59949709E-02 

70 

1.14405427E-07 

-1.96104461E-07 

1.02615469E-05 

6,28887408E-02 

71 

1.65745567E-07 

-8.73802977E-09 

-5.82660431E-09 

2.88428069E-02 

72 

1.92160659E-07 

2.30933170E-07 

-1.26321987E-05 

-4.70731332E-02 

73 

1.05031425E-07 

3.95379757E-07 

-2.15199960E-05 

-9.679524936-02 

74 

-4.90093406E“0e 

4.17721507E-07 

-2.26210193E-05 

-1.10635773E-01 

75 

-1.73529469E~07 

2.81843089E-07 

-1.48179466 E-05 

-8. 67100390 E-02 

76 

-2.34182913E-07 

9.58375632E-09 

6.39701603E-09 

-3. 15968803 E-02 

77 

-2.58710335E-07 

-3.20106295E-07 

1.74149936E-05 

6.96391649E-02 

78 

-1.33393533E-07 

-5,38503291E-07 

2.92702575E-05 

1,336613456-01 

79 

7.79210795E-08 

-5.50691973E-07 

3.03633453E^05 

1.47419054E-01 

80 

2.40603980E-07 

-3.70626319E-07 

1.95833395E-05 

1.08943509E-01 

81 

3.02289280E-07 

-8.52149244E-09 

-5.69499266E-09 

2.80636251E-02 

82 

3.16278412E-07 

4.04730815E-07 

-2.i8899554E-05 

-9.40429520E-02 

83 

1.52C77748E-07 

6.67749535E-07 

-3.62463455E-05 

-1.681091836-01 

84 

-1.1131X952E-07 

6.84191543E-07 

-3. 70558546E-05 

-1.78314582E-01 

85 

-3.03512382E-07 

4.42124012E-07 

-2.34938675E-05 

-1,230896866-01 

86 

-3.53028614E-07 

4.79054708E-09 

3.21738900E-09 

-1.57017051E-02 

87 

-3.45113853E-07 

-4.61421596E-07 

2.47843487E-05 

1.15363571E-01 

88 

-1.49747724E-07 

-7.43563747E-07 

4.02956732E-05 

1.90551988E-01 

89 

1. 449709 92E-07 

-7.46556218E-07 

4.04388090E-05 

1,924736136-01 

90 

3.44911857E~07 

-4.67981720E-07 

2.50584371E-05 

1,202622546-01 

n 

3,62973005E-07 

2.43286648E-09 

1,722267696-09 

-9.638418276-03 

92 

3.27532612E~07 

4.55146445E-07 

-2.42795296E-05 

-1,248701986-01 

93 

1,45090817E-07 

7.09547922E-07 

-3.85226372E-05 

-1.844317616-01 

94 

-1.11722331E-07 

6.97014670E-07 

-3.79640995E-05 

-1.764160326-01 

95 

-2.8I086829E-07 

4.28726627E-07 

-2,32819862E-05 

-1.08527967E-01 

96 

-2.99574280E-07 

-9.37O70924E-O9 

-5.52645553E-09 

-6.637656736-04 

97 

-2.70353608E-07 

-4.41484170E-07 

2.31230666E-05 

8.134838556-02 

98 

-7.46228255E-11 

-7.28756070E-07 

3.92369407E-05 

1.45920054E-01 

99 

5.21414819E-07 

-7.71134756E-07 

4.27008188E-05 

2.034724756-01 

100 

7.59850632E-07 

-4.59625262E-07 

2.88703579E-05 

1.97151230E-01 

lOI 

-4.25301149E-07 

0. 

0. 

0. 

NUMBER  OF  ITERATIONS^  28 


MINIMUM  TA*  5,36074232E^02  N=  5.70000000E<-01 


RUNNING  TIME  APPROXIMATELY  72  SECONDS 


APPENDIX  E 


DEVELOPMENT  OF  MORE  ACCURATE  RING  THEORY 

Near  the  completion  of  this  investigation,  Dr.  David  Bushnell  of  the  Lockheed 
Missiles  and  Space  Company  verified  the  results  of  some  of  the  present  calculations  by 
using  a  newly  developed  program  for  shells  of  revolution  (ref.  26).  In  doing  so,  he 
obtained  results  which  indicated  that  the  thermal-buckling  behavior  of  the  ring-stiffened 
shells  investigated  herein  could  be  influenced  by  the  type  of  ring  theory  employed  in  the 
analysis.  Preliminary  results  suggested  that  the  out-of-plane  ring  bending  stiffness 
during  buckling  was  especially  important. 

The  purpose  of  this  appendix  is  to  present  an  extension  of  the  present  analysis  to 
include  the  first-order  effects  of  out-of-plane  deformations  in  the  rings  during  buckling. 
The  ring  theory  is  developed  by  a  method  proposed  in  reference  27  to  obtain  a  consistent 
theory  for  sturdy  rings  attached  firmly  to  a  shell.  A  key  assumption  in  the  development 
is  that  compatibility  of  the  ring  and  shell  deformations  is  enforced  along  a  normal 
passing  through  the  centroid  of  the  ring.  Thus,  local  deformations  of  the  ring  elements 
are  ignored  and  the  ring  behaves  as  a  rigid  beam.  More  exact  models  for  rings  which 
enforce  compatibility  only  at  the  point  of  attachment  have  been  proposed.  (See,  for 
example,  ref.  28.)  However,  these  models  lead  to  extremely  complex,  unwieldy  equations. 

Ring  Displacements 

The  ring  center  of  shear  and  centroid  are  assumed  to  coincide.  The  ring  displace¬ 
ments  u|.,  vt ,  and  w*  are  specified  in  local  ring  coordinates  x,  y,  and  z,  which  are 

parallel  to  the  coordinates  of  figure  1  but  have  their  origin  at  the  ring  centroid.  The  ring 
displacements  u^ ,  v* ,  and  wt  are  expressed  (ref.  28)  as 

ut  =  Uc  +  Z)3j 

vt  =  Vc  +  z/32  -  xug -I  (El) 

wt  =  Wc  -  XiSj 
Vc 

with  =  -w^,  and  ^2  =  v  '''  where  the  subscript  c  denotes  the  ring  centroid. 

^  rCr* 
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The  ring  coordinates  y  and  z  are  related  to  the  shell  coordinates  y  and  z  by 


y=  i.-y 


z  =  z  -  z 


The  displacements  in  the  shell  in  the  vicinity  of  the  ring  are  assumed  to  be 


“sk  -  “  - 


''k  =  '  -  T,y  -  S 


wK.  =  w 


To  develop  a  ring  theory  for  sturdy  stiffeners,  compatibility  of  shell  and  ring  defor¬ 
mation  is  enforced  along  a  shell  normal  passing  through  the  centroid  of  the  ring.  Matching 
ring  and  shell  displacements  aloi^  x  =  0  /that  is,  ,  and  so  forthj  yields 


Uc  =  U  -  ZfW^x 


V,.  =  V  1  +-^  -  Zv.W  V  +  0 


=  w 


U*  =  U  -  ZW  5 
I  5-^ 


vt  =  V  1  +  -  .  zw^y  -  X  1  -  -  Ky  -  ZrW^jyU  0  - 


wt  =  W  +  XW  y 
r 


Strain  Energy 

The  strain  energy  in  the  N  interior  rings  on  the  shell  is  approximated  as 

"  I  Z  Iq  GrJr/sf  y  -  ErQ^TeJ.  6(x  -  jZ)dx  dy 
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t  ? 

t  t  1/  tV 


where  the  introduction  of  the  Dirac  delta  function  permits  the  replacement  of  x  by  x 
in  the  derivatives  appearing  in  equation  (E5). 

Buckling  Strain  Energy 

By  following  the  procedures  outlined  in  appendix  B,  the  total  potential  energy  of  the 
shell,  its  stiffening  elements,  and  the  external  loading  can  be  formed,  and  the  principle  of 
minimum  potential  energy  (611  =  0)  can  be  used  to  obtain  nonlinear  equilibrium  equations. 
The  equations  governing  prebuckling  and  buckling  can  be  derived  by  separating  displace¬ 
ments  into  axisymmetric  and  small  asymmetric  parts  (for  example,  u  =  u^  +  ug).  For 
the  more  accurate  ring  theory,  the  prebuckling  equations  obtained  are  the  same  as  those 
developed  as  equation  (1)  or  equation  (B21)  in  appendix  B.  Thus,  the  new  terms  of  the 
more  accurate  ring  theory  appear  only  in  the  budding  equations. 

In  order  to  make  a  comparison  of  the  present  theory  with  existing  ring  theories,  the 
strain  energy  during  buckling  can  be  written  as 


3=1 


where 


elx,  =  V 


rB  =  ^B,y  ^  r  "  ^^B,yy  “  ^  ^  "  RA^B,yy  ' 


xyy  +r^B+^B,x 


%  =  -^B.x 


(ElO) 


To  express  the  energy  as  a  function  of  the  buckling  displacement  variables  U,  V, 
and  W,  equations  (C13) 

TT  ny 

Un  =  U  cos  -r 
^  R 

Vd  =  V  sin  — 


Wt5  =  W  cos 
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are  employed.  The  integration  over  the  ring  area  can  be  accomplished  with  the  following 
definitions: 


^  Ep  dAp  —  EpAp  ^  EpX  dAp  —  0 

^  E-^z  dAr  —  ZjiEpA];*  ^  EpXZ  dAj*  = 

J  EfZ^dAr  =  Erioz  ^  Erx2dAr  = 


Xr  =  0 


dAf  -  Ei-Iqxz 


(Ell) 


where  the  subscript  o  is  used  to  denote  that  the  inertias  Iqz  and  Iqx2;  are  to  be  com^ 
puted  about  the  shell  reference  surface,  z  =  0. 

After  integration,  equation  (E8)  can  be  written  in  matrix  form  as 


in  which 


= r 

2  Jf 


Z^GZg  dx 


Zb  = 


(E12) 


(E13) 


Symmetric 


^12 

^13 

Gi4 

G22 

^23 

G24 

^33 

G34 

G44 

(E14) 
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with 


N 


Gii  =  ^  5(x  -  j 

j=l 


R 


N 


G12  =  2  -  i^)’7ErIo5z  h 

j=l 


N 


Gi3  =  2  6(x  -  il)vErloxz(^ 
j=l 


N 


Gi4  =  2 
j=l 


^Erl^ni  .  r,2z-rErl5(|'^ 
R 


N 


G22  =  J  6(x  -  j^) 

j=l 


/  \^  2  2 
ErArl^  1  +  2zrErAr  —  +  Erl,  ” 


.R 


T^OZ  , 

R*"  R 


N 


G23  =  2  -  j^) 

j=l 


ErAr  ^  +  ZrErArP^)  +  Erioz  4 

R^  V  R-^  /  R^ 


N 

G24  =  2  6(x  -  jZ)(ErIojj2  ^  ^ 

1=1  V  R  R 


]= 

N 


G33  =  J  6(x  -  jZ) 

i=i 


Ey.Ar  ^  j^2 

-^  +  2zrErArii^  +  ErIoz(- 
R' 


R' 


R 


N 


G34  =  ^  6(x  -  jZ) 
j=l 


N 


G44  =  6(x  -  jZ) 

j=l 


®1^^0XZ  J^3  “  ^r^^rloiczfRj 


E„I- 


R 


-  25r7?ErI-  +  G^J 


2  r  '  ^  X  R 


^r3 
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where 

=  1  —  ^Note  thnt  ^p^oz  ”  Epij*  +  z^Ej^Aj*^ 

In  this  form,  the  elements  of  the  ring  stiffness  matrix  G  can  be  compared  with  results 
obtained  in  other  buckling  studies  of  ring-stiffened  shells.  (See  refs.  28  to  30.) 

Equilibrium  Equations 

The  condition  for  vanishing  of  the  first  variation  of  the  ring  energy  =  oj  yields 

the  contribution  of  the  more  accurate  ring  theory  to  the  buckling  equilibrium  equations. 

If  the  terms  in  6U,  6V,  and  6W  are  collected,  the  set  of  buckling  equations  which  are 
the  counterpart  of  equations  (CIS)  with  a  more  accurate  ring  theory  can  be  written  as: 

AiiU"  +  ^Ai2  -  Giiju  +  AisV  -  G12V  +  (ai4  -  TaAi5)w"  +  (ai6  -  G14  -  TaA^Jw' 

+  ^Ajg  -  Gi3  -  TaAi9^W  +  AjioM'  =  0  (E16a) 

A21U'  -  G12U  +  A22V’'  +  ^A23  -  r2i)v  +  A24W''  +  ^A25  -  G24  -  TaA26]W’ 

+  (a27  -  r22  -  TaA28)w  =  0  (E16b) 

(A31  -  Gi4  -  TaA32)u'  +  r3iU  +  A33V"  -  G24V'  +  (A34  +  r32  -  TaA35)v 
+  (^36  -  r'33  -  Ta^37)w"  +  (A38  -  r34  -  TAA39  -  TiAgipjw- 
+  ^Agji  -  G34  -  TaA3J2)w  -  M"  =  0  (E16c) 

A41U' +  A42V  +  A43W"  +  (A44  -  TAA45jW' +  A46W  +  M  =  0  (E16d) 

where 

N  2 

^21  =  ^22-  Z  5(x  -  ji)ErAr/|) 
j=l  ^  ^ 

N 

^22  =  ^23  ■  Z 
3=1 
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^31  "  ^13  ■  ®ll4 


Equations  (E16)  were  cast  in  finite -difference  form  using  the  techniques  described 
in  appendix  C  and  reference  17;  the  results  were  programed  as  part  of  the  determinant 
plotting  solution  in  the  computer  program  entitled  BAMSOC  I  (see  appendix  D).  However, 
numerical  results  employing  the  more  accurate  ring  theory  are  not  available. 
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TABLE  I.-  STRUCTURAL  PROPERTIES  FOR  BUCKLING  CALCULATIONS 


(a)  SI  Units 


Type  of 

IS9 

Er 

Ir 

*- 

Zr 

Gr  Jj* 

Eg 

-^s 

^s 

*_ 

zs 

GsJs 

cylinder 

GN/m2 

mm 

GN/m2 

cm2 

cm4 

cm 

N-m2 

Ring- stiffened 

183 

1.43  X  10-5 

183 

0.626 

1.228 

-1,364 

0 

Stringer-stiffened 

72.4 

.32 

2.52 

72.4 

2,394 

3.663 

-2.344 

68.44 

.30 

0.9 

113 

1.587 

5.732 

9.999 

113 

-.319 

.1699 

-.864 

0 

(b)  U.S.  Customary  Units 


Type  of 

«sk 

Er 

^r 

Ir 

♦  _ 

Zr 

Gj.Jr 

Es 

IBS 

Zs 

GsJs 

cylinder 

psi 

Per  OF 

psi 

in2 

in4 

in. 

lb-in2 

psi 

■a 

Ring-stiffened 

26.5  X  10® 

0.3 

7.95  X  10“6 

26.5  X  106 

0.097 

0.0295 

-0.537 

0 

Stringer  -stiff  ened 

10.5 

.32 

14 

10.5  X  106 

0.371 

0.088 

-0.923 

23  800 

Ring-  and 

stringer -stiffened 

14.5 

.30 

5 

16.4 

.246 

.1377 

-1.992 

3477 

16.4 

.0495 

.004082 

-.342 

0 

*Based  on  reference  surface  at  skin  middle  surface. 


RING  DETAIL 

Figure  2,-  Geometry  of  ring-stiffened  cylindrical  shells  from  reference  11. 

Dimensions  are  in  inches  (cm). 
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Figure  3.-  Thermal  buckling  predictions  for  ring-stiffened  cylinders. 


Figure  4.-  Buckling  interaction  curve  for  ring-stiffened  cylinders  under 
combinations  of  axial,  compressive,  and  thermal  loadings. 


STRINGER  DETAIL 

Figure  5.-  Geometry  of  stringer-stiffened  cylinder  from  references  23  and  24. 

Dimensions  are  in  inches  (cm). 


A 


a 

(a)  Longitudinal  skin  stress. 

Figure  6.-  Stress  distribution  in  clamped  stringer-stiffened  cylinder  under 
combinations  of  axial  compression  and  thermal  loadings.  (Negative  sign 
on  stresses  denotes  compression.) 
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Figure  6.-  Concluded. 


X  LOADINGS  OF  FIGURE  6 


300  (762) 


-40 


a 

(a)  Skin  stresses. 


Figure  10.-  Stress  distribution  in  simply  supported  ring-  and  stringer-stiffened 
cylinder  under  a  combination  of  axial  compression  and  temperature 
=  1650  Ib/in.,  0.289  MN/m;  Tgk  =  360°  F,  200  k). 
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Figure  10.-  Concluded. 


- AXIAL  LOAD  APPLIED  AT  SKIN  MID-PLANE 

.oo4p - axial  load  applied  at  skin-stringer  neutral  surface 
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combinations  of  axial  compressive,  thermal,  and  pressure  loadings. 


NASA-Langley,  1971  -  32  L-7021 


Figure  12.-  Radial  bucklins:  deformation  in  thermally  loaded  ring;-  and  strinerer-stif 
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